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Emptiness formation probability of the six-vertex model and 
the sixth Painleve equation 

A. V. Kitaev and A. G. Pronko 


Abstract. We show that the emptiness formation probability of the six- 
vertex model with domain wall boundary conditions at its free-fermion point is 
a r-function of the sixth Painleve equation. Using this fact we derive asymp¬ 
totics of the emptiness formation probability in the thermodynamic limit. 
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1. Introduction 


One of the most intriguing facts about correlation functions of solvable models 
of statistical mechanics is that in many cases they can be described in terms of 
of the Painleve equations or their generalizations. Famous examples are provided 
by correlation functions of the two-dimensional Ising model, related to the third 
Painleve equation mm, and correlation functions of an impenetrable Bose gas, 
related to the fifth Painleve equation m- At the same time this relationship carries 
mainly academic interest, whilst asymptotic analysis of the correlation functions is 
usually performed via the associate integrable structures rather than with the help 
of these ordinary differential equations (ODEs). 

In this paper, we consider the six-vertex model with the domain wall bound¬ 
ary conditions and discuss a particular correlation function, called the emptiness 
formation probability (EFP). We show that, for the model with the Boltzmann 
weights satisfying the free-fermion condition, this correlation function appears to 
be a r-function of the sixth Painleve equation (P6). Using this connection we derive 
asymptotic expansions of the EFP in the thermodynamic limit. 

It is important to stress that here we consider asymptotics of P6 for the case 
where the coefficients of the equation are large while its argument is a finite param¬ 
eter. Such type of asymptotics is not thoroughly studied for P6 yet. Therefore to 
cope with this problem, we initially intended to use one of the known asymptotic 
methods for the Painleve equations, namely, the isomonodromy deformation tech¬ 
nique by Jimbo-Its-Novokshenov or Deift-Zhou [7] asymptotic analysis of the 
corresponding Riemann-Hilbert problem. However, we were quite surprised that 
we have been able to construct asymptotics (both the leading terms and correc¬ 
tions) without any sophisticated techniques—^just with the help of the cr-form of 
P6. Indeed, finding asymptotics, especially the correction terms, of solutions of 
ODEs by substituting corresponding asymptotic expansions into the ODE is the 
standard asymptotic technique in the case where the asymptotics with respect to 
the argument of the corresponding ODE is constructed. On this way, one usually 
finds a recurrence relation for determination of the coefficients. Our case is differ¬ 
ent: we arrive at the recurrence relation for the first derivatives of the correction 
terms so that to find the explicit formulas we need the initial data. It turn out 
that for correction terms of arbitrary order these initial data can be obtained from 
representations of the EFP at the critical points of P6. It seems that the same 
technique may appear useful for other analogous problems. 

1.1. The model. We recall that the six-vertex model (also known as the ice- 
type model) is a statistical mechanics model defined on a square lattice. The local 
states of the model are arrows placed on the edges of the lattice. The admissible 
conhgurations of the model are those in which there are equal numbers of incoming 
and outgoing arrows in each lattice vertex. The condition of local conservation 
of the number of incoming and outgoing arrows is known as the ice-rule, and it 
selects six possible arrow configurations around the vertex, see Fig.[T] For periodic 
boundary conditions the model was solved by Lieb [BHin] and Sutherland [H] (for 
a review, see, e.g., [1^1. 

Domain wall boundary conditions (DWBC) are defined as follows. Consider a 
finite square lattice formed by intersection of N horizontal and N vertical lines (the 
so-called N x N lattice). For such a lattice one may require that the local states on 
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Figure 1 . The six types of arrow configurations around a vertex 
allowed in the six-vertex model and their Boltzmann weights. 


s + q r 



Figure 2. An N x N lattice {N = r + s + q) with domain wall 
boundary conditions, with the additional condition that the ver¬ 
tices of an s X (s -I- q) lattice at the top left corner are all of type 
2—the probability of this configuration is Fr^s,q (shown the case 
r = 5, s = 3, and q = 1). 


external edges are fixed in a special way consistent with the ice-rule, namely, using 
the description of the local states in terms of arrows on edges, that all arrows on 
the external vertical edges (i.e., at the top and bottom boundaries) are incoming, 
while those on the external horizontal edges (i.e., on the left and right boundaries) 
are outgoing, see Fig. [2] The six-vertex model on the N x N lattice with such fixed 
boundary conditions is called the six-vertex model with DWBC; it was considered 
for the first time by Korepin m- 

To define the partition function of this model, let us denote by the set of 
all arrow configurations of the N x N lattice with DWBC and with the ice-rule 
obeyed in each lattice vertex. Assigning to the ith. arrow configuration around a 
vertex the Boltzmann weight Wi, see Fig. [U the partition function Zn is defined as 

= E n 

Cg^n 4—1,...,6 

Here ni(C), i = 1,..., 6, is the number of vertices of type i in the configuration 
C, = N'^. The partition function was computed, for a more general 

inhomogeneous (with position-dependent weights) model, in terms of an A^ x A^ 
determinant by Izergin |14j : for the homogeneous model the determinant has the 
Hankel structure m- The free energy per site was derived for various regimes 
and by various methods by Korepin and Zinn-Justin [mill]. A detailed analysis 
of the asymptotic expansion of the partition function in the thermodynamic limit, 
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using the connection of Hankel determinants and matrix models with the Riemann- 
Hilbert problem [18] . was given in paper |19j and a series of papers by Bleher and 
collaborators [5DH23]. 


1.2. Emptiness formation probability. To define the EFP, let us set N = 
r + s + q, where r, s, and q are nonnegative integers, and consider configurations of 
six-vertex model on the N x N lattice with DWBC in which the vertices belonging 
to the s X (s -b g) rectangle at the top left corner of the lattice are all of type 
2, see Figs. [Hand [5] We denote the set of such configurations as ^r,s,q', obviously, 
rir-.s.g C Oat. We denote the EFP by Fr^s,q and dehne it as a probability of observing 
a configuration belonging to flr,s,q for a given r, s, q, 

Fr,s,q=Z^^ E n 

CGOr,s,g i—1,...,6 

In what follows we always assume that s < since otherwise Fr^s,q is jiist zero 
(for s > r the set ^lr,s,q is empty). A closed expression for the EFP in the terms 
of a multiple (s-fold) integral was obtained by Colomo and the second author in 
[25] . As it was demonstrated in [26], this result make it possible to address phase 
separation phenomena in the model, and in one of its special cases provides the 
limit shape of large altering-sign matrices m- 

In this paper we study the EFP of the six-vertex model with DWBC in the case 
of the weights obeying the so-called free-fermion condition, wiW 2 + 'W 3 W 4 = w^wq. 
Namely, we choose the weights in the form 

wi = W 2 = — a, W 3 = 104 = Va, W 5 = wq = 1 , ( 1 - 1 ) 

where a G (0,1) is a parameter. We note that this case is related to enumeration 
of domino tilings, where the parameter a plays the role of “bias” [251130] . In the 
main text we prove the following result. 


Theorem l.I. Define 


I d 1 r (rFq + sf 

a = a{a - 1)— logI<r,s,q -- a 

da 4 


(r + g + s)q + 2 rs 


( 1 . 2 ) 


then 


a^{a — -I- {(1 — 2a){a')‘^ + 2aa + i'ii'2t^3i^4:}^ 

= + ^ 1 ) + ^ 2 ) + ^ 3 ) + ^' 4 ) > cr' = Ecr, (1.3) 

where 

r-b g-b s r — q — s r + q — s 

U,=U3 = - ^2 =-= . 

Remark 1.1. Note that Fr^s,q is a polynomial in a, see for details represen¬ 
tation (12.21) obtained in Sect. I^.IL Thus, m implies that a = (7{a) is a rational 
solution of m. The corresponding solution of the canonical form of P6 is also 
rational. This solution defines isomonodromy deformations of the associated linear 
Fuchsian ODE with the monodromy matrices ±/ (see Sect. 15'.31 for details). 


Remark 1.2. The EFP can be expressed as 

= (1.4) 

where quantity Cr,s,q is independent of a and t = ria) is the Jimbo-Miwa r- 
function of P 6 [4], see also Sect. fOl The pre-factor of the t- function in (d is 
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a matter of the definition (normalization) of the r-function, therefore Fr^s,q is, in 
fact, the T-function of P6. 

We recall some basic fact about the theory of P6 in Sect. 13.11 Our starting 
point in proving the theorem is the formula for the EFP expressing it as a Fredholm 
determinant of some linear integral operator of the so-called integrable type. This 
representation, among some others, was derived from the general multiple integral 
representation by specifying it to the case of the weights (EB in |3I]. We review 
some of these results in Sect. 12.1] The theorem is proven in Sects. H] and [3] 

As a comment to this result, we mention that appearance of a Painleve equation 
in the context of the six-vertex model with DWBC was already noticed in [32] where 
it was shown that the partition function in the ferroelectric regime serves as a r- 
function of the fifth Painleve equation. In turn, P6 appeared in the study of the 
FFP at zero-temperature for the Heisenberg XXO spin chain [33] , which can be seen 
as the six-vertex model at its free-fermion point with periodic boundary conditions. 

In fact, the relation of Fr^s,q with P6 reported in Theorem 1 1.1 1 can be extracted 
from the works by Johansson, Baik and Rains, and Forrester and Witte [33H33] (see 
a discussion at the end of Sect. 1311) . Here we present a different proof based on 
Fredholm determinants of integrable linear integral operators, developed in [31137]. 
This method is also important for asymptotic study of Fr^s,q as it can be used to 
provide initial conditions in our proof of Theorem ll. 31 below. The connection of the 
FFP with P6 allows us to study behavior of the FFP in the thermodynamic limit, 
in which the integers r, s, and q are large, with their ratios fixed; the variable a, 
defining the weights, is considered as a parameter. An important special case of 
such a limit is where the region of the vertices of type 2 at the top left corner has a 
macroscopically square shape, that corresponds to the ratios q/s and q/r vanishing 
in the limit. In this paper we compute the asymptotic expansion for Fr^s,q in the 
case of g = 0. 


1.3. Main result: asymptotic expansions. Define 




v,u G (0,1). 


As discussed in Sect. oi the FFP has different asymptotic expansions for u < u < 1 
(referred to as the FFP in the disordered regime) and for 0 < u < u (referred to as 
the FFP in the ordered regime). Our results about the FFP in the thermodynamic 
limit are formulated as follows. 


Theorem 1.2. For s/r = v G {u, 1), as s,r ^ oo. 


nil (1 — u^)u^ 1 1 — , 

log F,,,o = -».»- 1“8 »+ 8 >«« - T2 ^ 


+ 0 -2kS -f 0(s ^), (1.5) 


fc=l 


where 


if = log-- 
u 


(1 — vY 1 — u (1 -f vY 1 + 


■log 


log- 


( 1 . 6 ) 


1 — u 2^2 1 -I- u ’ 

C'(—1) = —0.1654211437... is the derivative of the Riemann zeta-function, and 
a 2 k = a 2 k{u,v) are rational functions ofu^ andv^ with poles atv^ = 1 andv^ = v?. 
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In particular, 

u'^{\ — v"^){2v‘^ + — u^) {1 + v'^){v'^ — {1 — v"^Y) 1 

^ 64(^2 - u 2)3 120(1-i;2)2 

+ 140i;®u^ + 105?;®u® - 160i;®w'^ - 4v*u^ + - 6v^u^ + + u^°] 

v^{v^ - 4^4 + 5z;2 - 10) , 31 

504(1-?;2)4 16128’ 

The estimate 0(s“^"“^) of the error term in (II.5p is uniform with respect to the 
parameter v on any compact subset of the interval {u, 1). 


Theorem 1.3. For s/r = v G (0,m), as s,r ^ oo, 


log (1 - Fr,s,o) = -xs - log s + log 


^(1 — M^)2 


32tt^/u[{u — z;)(l — uv)]^/'^ 

n 

+ ^6fcS-'= + 0(s-’^-i), (1.7) 


k=l 


where 


4 , / VI — vu + \/u{u — v) 

* = — 77=^ — 


- 4 log 


\/u{l — vu) + y/u — ‘ 
y/{l- U^)v 


and bk = bk{u,v). In particular, 

9?;^(1 + u'^) — 36vu{l + f^)(l + u^) — (2 — 142ti2 + 2v^)u‘^ 


bi = -- 
&2 = 


48Vw[(u — z))(l — uv)]^/'^ 

■ [3t>^(u® + 1) — 8vu{v‘^ + l)(u® + 1) 


( 1 . 8 ) 


64m(m — ti)^(l — vu)^ 

+ 4M^(10ti^ + v‘^ + 10)(u'^ + 1) ~ 120vu^{v‘^ + l)(u^ + 1) 

- (16i;^ - 370?;2 + 16)u^]. 

The estimate 0(s“"“^) of the error term in (11.71) is uniform with respect to the 
parameter V on any compact subset of the interval (0,u). 


Remark 1.3. Since s and r are integers v is a rational number, v = so/ro, 
where 0 < sq < ro are co-prime natural numbers. Asymptotic expansions pre¬ 
sented in Theorems \1.S\ a,nd \l.S\ are understood in the natural way: put s = Son 
and r = ron, where n is a positive integer and n —>■ oo. As a further comment, 
we note that one can modify the problem and ask about asymptotics for irrational 
values of V. Our results also deliver the answer to this question. Consider rational 
approximation, v = Sn/rn + Sn, where Sn andrn are co-prime natural numbers, and 
Sn ^ 0 is the corresponding error. With the help of the Farey series [3S] we can find 
such approximations that Sn = where |k„| < 1. For the irrational v, Theo- 

rems \l.‘A and \l.A remain valid but with v substituted by v — Sn- In particular, for the 
optimal rational approximations (obtained with the help of the continued fractions 
or the Farey series) it is convenient first to rewrite the asymptotics in terms of the 
large parameter r instead of s by using the relation s = vr, then substitute v with 
V — Knfrf^, and, finally, re-expand asymptotics with respect to r = Vn. 
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Remark 1.4. Existence of the asymptotic expansions dm and (Ell are proved 
in Seet. \4.‘d\ and Sect. respectively, and the eoefficients a 2 k’s and bk’s can be 
obtained by recurrent procedures via substitution of these expansions into the a- 
form of P6, see Theorem 11.11 The coefficients are derived successively with the 
linear differential equations of the first order whose initial data are given in these 
sections. 

The asymptotic formulas (USD and EZl confirms and extends the result of 
[30] , where only the leading 0(s^)-term of the asymptotic expansion of log^V.s.o 
was obtained. Formula (11.81) for y reproduces the result obtained by Johansson [34] . 
see, for details, Sect. oi 

As a comment to the asymptotic expansion of the EFP, we mention that the 
value V = u is the critical value at which the third-order phase transition of the 
leading term takes place. This phase transition has a combinatorial interpretation 
in the context the of random domino tilings on Aztec diamonds with a cut-off corner 
m- Recently, a possible third-order phase transition in the thermal correlations 
of XXO spin chain, although for a peculiar choice of the parameters, was pointed 
out in [39] . Similar phenomena can be seen when the EFPs of quantum spin chains 
and vertex models are treated by conformal field theory methods [40] and when the 
vector Chern-Simons theory is considered on x [41]. A detailed discussion of 
the third order phase transition phenomena from the random matrix models point- 
of-view can be found in [42]. We discuss aspects of the third-order phase transition 
related to the EPP in terms of the associated Painleve equations in Conclusion. 


2. EFP and integrable structures 

The aim of this section is to present an overview of the results which lead to 
formulation of the EFP as a special solution of an integrable system (in our case 
P6). We start with recalling presentation of Er.s.g as a Fredholm determinant of a 
linear integral operator, next we discuss the properties of the corresponding resol¬ 
vent kernel, and finally we show how to convert the associated system of integral 
equations to Fuchs (a Lax-type) pair for P6. 


2.1. EFP as a Eredholm determinant. We mention here three represen¬ 
tations for the EFP, valid for the six-vertex model with the weights m-- in terms 
of a multiple integral, in terms of Hankel determinant, and in terms a Fredholm 
determinant of a linear integral operator acting on a closed contour in the complex 
plane. Others can be found in m where interrelations between these representa¬ 
tions are proved. In our study of the EFP we mostly rely on the last form, in terms 
of the Fredholm determinant. 

In [25] , using the method of Yang-Baxter algebra (the algebra of the quantum 
monodromy matrix) [43] and the methods of orthogonal polynomials theory, the 
EFP was evaluated in the form of a multiple integral. In the case of the weights 
CH) the result can be formulated as follows. 


Proposition 2.1. The EFP of the six-vertex model with DWBC with the 
weights EB is represented in the form 


jr — 

^ r,s,q — 


(_l)s(s-l-l)/2 

.s! 


Co 



{azj -I- 1 — d®z 

z'jizj-iy (27ri)«’ 

( 2 . 1 ) 
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where Cq is a simple, closed, counter-clockwise oriented contour around the point 
z = 0 of the complex plane, and lying in its small vicinity. 

The representation (12.11) was used in the study of the so-called arctic curve 
of the six-vertex model with DWBC, i.e., the curve which describes the spacial 
separation of order and disorder in the thermodynamic limit [261144] . However, the 
integral can hardly be studied by the random matrix model methods, since 
the problem of finding of an equilibrium measure cannot be solved explicitly (see 
discussion in [44]). This fact stimulated to look for another, simple representations. 
In particular, evaluating the integrals it can be shown that the EFP is given in 
terms of a Hankel determinant. 


Proposition 2.2. The EFP admits the representation: 

' r — 1 


Er.s.o — 


iq\r 


(1 - 


UlJoiQ + ky.kl 


q,s(s-1)/2 


det 

l<j,k<s 


E 

rn—O 




( 2 . 2 ) 


Using a peculiar structure of the matrix in (12.2[) . the finite-size determinant here 
can converted into Fredholm determinants of integral operators of various types. 
Indeed, the matrix in (12.21) can be represented as a difference of two matrices, one 
involving the summation over m G Z>o and second over m G Z>j.; the first matrix 
is related with the ensemble of Meixner polynomials and hence it can be explicitly 
inverted, while the second one gives rise to a Fredholm determinant structure of the 
finite-size determinant. This finite-size determinant can be further transformed by 
keeping the Fredholm structure of the determinant. In particular, it can be written 
as a Fredholm determinant of a linear integral operator acting on a closed contour 
in the complex plane. Specifically, here we deal with the integral operator K acting 
on trial functions by the formula 

(k/) (A)=£ 

The function K(x,y) is the kernel of the operator, and Cq is the contour given as 
a circle of small radius around the point A = 0, with counter-clockwise orientation, 
i.e., the same contour appeared above in the multiple integral representation for 
the FFP. Besides this contour, the functions defining the kernel involve integrals 
over the contour Coo, which is a circle of large radius, counter-clockwise oriented 
around the origin (the circle of a small radius around the point A = oo, clockwise 
oriented around it). 

We also recall that the Fredholm determinant of a linear integral operator is 
defined as 


Det (l - A') = 1 + f; ^ *t^^ [A(A., a,)] d"A. 

n—1 

Using that Det(l — K) = exp{Trlog(I — K)}, and defining the function log(l — K) 
by its power series expansion in powers of K, we can also define the Fredholm 
determinant of the integral operator K by the formula 

Det(l-i:) =exp|-£i J K(Ai,A 2 )---K(A„,Ai)d"A|. (2.3) 

We now formulate the Fredholm determinant representation for the EFP, derived 


in 






Proposition 2.3. The EFP admits representation in terms of the Fredholm 
determinant 

(2.4) 


Fr,s,q = Det ( 1 - 


(l - A-) 


where K = Kr,s,q is the integral operator acting on the contour Cq and possessing 
the kernel 

e+(A)e_(/i) - e_{A)e+(/^) 


KiX,fi) = 


X, 1-1 G Co 


(2.5) 


27ri(A — fr) ’ 

where the functions e±{X), which depend on the integers r, s, and q, i.e., e±(A) = 
^ih,r,s,g(A), are 

(A - ay/^ 


3_(A)= 


and 


where the function E(X) is 

Fix) = j 


(A - 1)C+<?)/2 as/2 
e+(A) = e_(A)£;(A), 

{v - ly+tv’^ die 
Coo ~ cxy{n — A) 27ri 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


and Co and Coo ci,re counter-clockwise oriented circles of small and large radii, 
respectively. 


The details of the proof can be found in m- Here we just mention the idea of 
the proof: upon change of the variable v i—^ l/p, which replaces the contour Coo by 
the contour Co, the trace of the nth power of the operator K reads 


Tr = (-!)”« (b 

Co 


where Ap+i := Ai. Since all A^’s and Vi are integrated over the contours of small 
radii, the factors (1 — XiVi)~^ and (1 — Xi+iVi)~^ can be expanded in the Taylor 
series. Re-ordering the summations and integrations here, various Fredholm-type 
determinant representations can be obtained, in particular, in terms of a finite-size 
matrix, which leads to the Hankel determinant representation (12.2p . 

It is necessary to mention that the representation (12.21) . up to a redefinition of 
the discrete parameters and when rewritten in its random matrix model integral 
form (with a discrete measure), appears to be coinciding with certain quantity 
discussed in a different context by Johansson [33] (see Proposition 1.3 therein). In 
turn, Baik and Rains |351I45] showed that this quantity can be represented as certain 
average over circular unitary ensemble. Using the Okamoto theory [46] . Forrester 
and Witte [36] next showed this average serves as r-function corresponding to 
a classical solution of P6. In what follows we show that the this result, which 
can be summarized in the form of Theorem 11.11 can be established directly from 
the Fredholm determinant representation (12.411 . using the notion of the integrable 
integral operators introduced by Its, Izergin, Korepin, and Slavnov [37]. 
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- iy.Y+<l 


{a - Ai)’^(l - i^i) 

Cot\ (l-A,)^+«(l-ap,)^ 


d"i/dPA 


Af+i<+''(l - A,p,)(l - A.+iPz) (2^i) 


2n 








2.2. The resolvent operator. Consider Fredholm integral equations 
/±(A)-/ K{X,^i)f±{n)dfj, = e±{X), 

JCo 

which can be written in operator form as 


(l-iF)/±] (A)=e±(A). 


(2.9) 


For brevity we omit dependence of all functions on the parameter a and discrete 
parameters r, s, and q. We consider solutions of (12.9p belonging to the class of 
analytic functions in A-plane which can be presented in a vicinity of the point 
A = 0 in the form 

/±('^) = —+ h2(A)e_(A), 
e_(A) 

where hi (A) and h 2 (A) are analytic functions at the origin. For given values of r, 
s, and q there exists a finite number of points a = ai,..., ap, where p = r{s + q), 
such that dimKer(l — K) = 0 for a G C \ {ai,... ,ap}, and dimKer(l — K) = 1 
at a G {«!,..., ap}. It is important to note that ai,... ,ap ^ (0,1) that is the 
operator 1 — IF is invertible for the case of our interest, a G (0,1). Therefore, one 
can define the resolvent R of the operator K, 


(i-itr) (i + -R) = 1 - 


( 2 . 10 ) 


As we show below, the kernel the operator R can be written in terms of the functions 
exactly in the same way as the kernel of the operator K is given in terms of 
the functions e±(A) [37] . 

It is easy to see that (12.9|) can be rewritten in the form 

X(A)e(A)=/(A), 


where 


and 


?(A) = 


(e+W\ 

V-{X)) 


/(A) = 


(UW\ 

[f-W ’ 


XiX) = 

where the functions Hmn{X) are 

Rnin{X) = (p 

JC 


l-i7_+(A) H++{X) 

-i7__(A) l + i7+_(A)r 


Cm(M)/n(M) d/r ^ 

-r —, TO,n={ + ,-}. 

Co h - A 27ri 


( 2 . 11 ) 


( 2 . 12 ) 


It turns out that the matrix Ar(A) can be easily inverted, due to the following 
remarkable property. 


Proposition 2.4. The determinant of the matrix (12.111) is equal to one, 
detX(A) = 1, A G C. 


Proof. To prove the result we use the same trick as in [33], Chap. XIV, 
App. C, but we will rely only the known structure of the kernel of the operator K, 
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instead of that for the resolvent. Note that the functions _(A) and H _^(A) can 

be written in the form 


and 


H+-iX) = 


H-+{X) = 


(l-x)/+ (Ai)/-(/r) 


Co 


27ri 


/i — A 


27ri’ 


I Co V -X 

respectively. Subtracting the second expression from the first one and using the 
fact that the kernel is symmetric, v) = K{y^ /r), we have 


- H-+(A) = £ * £ /-(rtK(/., 0/+(0 (;;:4 - 

Since 

1 1 fj, — 1/ 

v — X ^ — X (/r—A )(ja —A) 

and 

- v)K{^l, v) = e+(/r)e_(zA) - e-{^l)e+{v), 
it is fairly easy to see that in fact we have the relation 


7J+_(A) - H_+{X) = iI„+(A)iJ+_(A) - i?__(A)iJ++(A). 


Obviously, the last relation implies that detX(A) = 1. 


□ 


Corollary 2.1. The inverse matrix has the form 


x-\x) = 


fl + H+_{X) 

1 H..{X) 


-H++{X) \ 
1-H.+ {X)J- 


(2.13) 


Corollary 2.2. The functions e±(A) solve the following integral eguations 


f±{X) = e±(A) + (f R{X, At)e±(^) d^ 

JCo 


where 


R{X,ti) = 


/+(A)/_(Ai)-/_(A)/+(A.) 


27ri(A — fx) 

Proof. Using (12.121) it is straightforward to rewrite the formula 

X-i(A)/(A)=e(A) 

in the form of the integral equations (I2.14|) with the kernel (12.1511 . 


(2.14) 

(2.15) 


□ 


Note that (12.141) can be written in the following operator form 


/±(A)= [(l + i?)e±] (A). 


(2.16) 


Comparing (j2.16|l with p.9ll . we conclude that the operator R is the resolvent of 
the operator K, defined by the formula ()2.10l) . In other words, we have just showed 
that if the operator K possesses the kernel (12.5L then the functions /±(A), defined 
as the solutions of the equations (EH), determine the kernel of the resolvent by 
formula (12.151) . 

In a similar manner one can calculate action of the operator (1 — K)~^ = 1 + ii 
on various functions. In what follows we often use the following result. 
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then 


Lemma 2.1. Let the vector g(\]v) to he 

(i + i?)5(-;^)] (A) 
’roof, a direct calculation gives u 
(l + i?) 9±M] 


1 


X — V 

e±(A) _ 

X — V ^ X — V 

/±(A) 

\ — u 

/±(A) 

A — 


/\ \ 

9±{A,iy) = 

X — u 

(2.17) 

= y^X-1(p)/(A). 

A — ly 

(2.18) 

d/i 

3 T-v 



(f R{X, /i)e± (^) (^1 + ^ dfi 
Jco \ L-^J 


I Co 

X-u -R(A,/r)(A -/r)e±(/i)d/r 


X- 


X- 


(2.19) 


where relations (I2.16P and (I2.15P have been used. Comparing the result in ()2.19l) 
with (I2.13L we arrive at (I2.18p . □ 


2.3. The Fuchs pair. One of the main properties of the functions /±(A) is 
that that obey linear first-order differential equations in the variables a and A. 
These equations appear to be a vector form of the famous Fuchs pair for P6. 

Note that there exist also difference equations with respect to the discrete 
parameters r, s, and q, shifting them by ±1. These equations, together with the 
one with respect to a, constitute the set of the Lax-type equations of our problem. 
They can be viewed as the so-called Schlesinger transformations for the Fuchs pair 
given below, or derived in a straightforward way using the method presented in the 
proof. 

Proposition 2.5. The functions f±{X) is the solution of the system of linear 
differential equations 

A/(a) = A(A)/(A), A/(a) = B(A)/(A), 

with 

= X + + B(A) = -^, (2.20) 

A A—lA — a A — a 

where the matrices Aq, Ai, and A^^ are given by 

A, = ^X{,y)S{^)X-\iy) (i/ = 0,l,a). (2.21) 

Here 

6>o = s, ei=r + q, 9^ = -r, 
and the matrix X(A) is given in (I2.11[) and the matrix S{X) is 


( 2 . 22 ) 

































Proof. Here we consider a proof based on the straightforward differentiation 
of the functions /±(A); a simpler proof, based on analytical properties of the matrix 
X(A), is given in Sect. 13.11 

We first consider the derivatives of functions /±(A) with respect to the param¬ 
eter a. We begin with calculating the derivatives of functions e±(A). From (12.6L 
for the function e_(A) we have 



r 

2(A-a) 


e-(A). 


(2.24) 


To obtain the similar relation for the function e+(A), let us consider first the function 
if (A), see (12.71) . for which we have 


da^^ ^ "" Jc^ {v - aY+^iy - X) 2Tri' 


Using 


1 


1 


1 


{v — a){v — A) X — a \ ly — X 
we can rewrite the relation above as 

d 


-EiX) = - - 

da A — a 


(EiX) - E{a)) 


and therefore 

Ae+(A) = ^_e+(A) - ^e_(A), (2.25) 

da 2(A — a) A — a 

where we have also used (12.241) . 

To obtain derivatives of the functions f±{X) with respect to a, we differentiate 
the defining relation (12.91) . that gives 


d 

da 


/±(A) - 



da 


/±(m) 


Jco da 


= ;pe±(A). 

da 


(2.26) 


Using (12.241) and (I2.25L the derivative of K{X, Y) can be computed, with the result 


d e+{X)e-{ii)+e-{X)e+{^i)-2E{a)e-{X)e-{p) 

- 4,i(A-o)(^-») -■ 

Thus, moving the second term in the left-hand side of (12.261) to the right, we can 
rewrite (12.261) as follows: 




(A) = —e±(A) 

da 

■ [ii_±(a)e+(A) -I- ii+±(a)e_(A) — 2if(a)ii_±(a)e_(A)]. 


2(A-a) 

In the case of the plus-sign subscript, using (I2.25L we thus obtain 




(A) = (l-i7_+(a)) 


^e+(A) 

2(A-a) 


re(A) 


- [H++{a) + 2E{a){l - H_+{a))] 


(2.28) 
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Similarly, in the case of the minus-sign subscript, using (12.241) . we get 




(A)=-7J__(a) 


re+{X) 

2{X-a) 


- [1 + H+.{a) + 2E{a)H._ia)] (2-29) 

Now the desired result can be obtained by acting with the operator (1 — K)~^ = 
1 -I- i? on these two relations. 

Comparing (I2.28P and (12.291) with (12.111) , one can note that these two relations 
can also be written as 

d 




iX) = -Xia)Sia)g{X;a), 


where the matrix S'(A) is defined in (|2.23p . Specifying z/ = a in (12.181) . we obtain 

^/(A) = ^^j^X{a)S{a)X-\a) /(A). (2.30) 

Clearly, the matrix standing in the right-hand in (12.301) side is exactly the matrix 
B{X) appearing in (12.201) . 

Let us now find the derivatives of the functions /±(A) with respect to the 
variable A. We first obtain those of the functions e±(A). It is straightforward for 
e_(A), 

r 


^"-(^^-“2 Ia + A ^" A ^ 


e_(A). 


(2.31) 


To obtain that of e+(A), we note that the derivative of the function E{X), after 
integration by parts, reads 

iy%iy-lY+i diy 


3 a^(^' = 


Coo ~ aY{v — A)^ 27ri 

d i/®(z/-!)’'+«A 1 dv 


\dz/ Y — olY J v — X 27ri 

= Ue{X) - Em + ^(i?(A) - E{1)) - ^{E{X) - E{a)). 
A A — 1 A — Of 


Hence, 


d ^ f s r + q r 

dA®+^ ^ “ 2 I A A- 1 “ A-a 


e+(A) 


sE{0) (r -I- q)E{l) rE{a) 


X 


A- 1 


A — a 


e_(A), (2.32) 


where (I2.31|) have been used. 

Now we are prepared for calculations with the functions /±(A). Differentiate 
the defining relation (12.9p to obtain 




:jT + ^(A>A^) /±(M)d/r- K{X,n)-^f±{q)dq 


Co L\dA dfj. 


= —e±(A). (2.33) 
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Here the second term can be expressed in terms of the functions e±(A). Indeed, 
since 

(A. d \ 1 

\dX dfi) X — fj, ’ 

differentiation of the kernel affects only the numerator, 

M e'+{X)e-{fi) - e'_{X)e+{n)+e+{X)e'_{fj,) - e-{X)e'+{ii) 

U + 5?j ‘ -' 

where prime denotes the derivative with respect to the argument of the function. 
Using (12.311) and (12.321) . we obtain 


_ 1 
“ ~2 


s r + q r 

A/r (A — l)(/r — 1) (A — a){^ — a) 


+ 


sE{0) {r + q)E{l) 


rE{a) 


e+{X)e-{^i) + e-(A)e+(/r) 
27ri 

e-{X)e-{fi) 


Xfj, (A —l)(/r —1) {X — a){p — a) 


27ri 


Thus, moving the second term in the left-hand side of (12.331) to the right, we arrive 
at the following relation: 


(l - k) /i] (A) 




2(A-1) 


^ [H_±(0)e+(A) + (H+±(0) - 2E{0)H.±{0))e.{X)] 
[iJ_±(l)e+(A) -f (i?+±(l) - 2U(l)iJ_±(l))e_(A)] 


+ 2{x'’-a) + (^+±(«) - m(^)H-±{a))e.{X )]. (2.34) 

Using ()2.31l) and (I2.32p and switching to the vector form, one can see that ()2.34l) 
is just the relation 


(l - k) f'] (A) 


= I XiO)SiO)g{X-, 0) + ^ X(1)5(1)5(A; 1) 


^X{a)S{a)g{X\a), 


where functions g±(A; v) are defined in (12.171) . 

Finally, acting with the operator (1 — K)~^ = 1 -|- i? on the last relation by 
making use of p.l9p we obtain 


i/(A) = {^X(0)S(0)X-(0) + ^X(I)S(1)X-(I) 

- e(A). (2.35) 

Clearly, the matrix appearing in the right-hand side of (12.351) coincides with the 
matrix 7l(A) given in (12.201) . □ 
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3. EFP and P6 


The main goal of this section is to prove Theorem ll.il We begin with recalling 
some basic facts of the theory of P6. Next, we discuss analytic properties of the 
matrix X(A) in terms of the related objects of the theory of integrable systems. 
Finally, we use all these ingredients to show that the Fredholm determinant of the 
operator K serves as a r-function of P6. 


3.1. The sixth Painleve equation. The canonical form of P6 (cf. [47]) 
reads. 


da 2 


1 


1 


y-l y-a 

y{y- ^){y-o‘) 
— 1 ) 


dy 

da 


- + ) — 

a a—1 y — a/da 

Ma - 1) 

(y - a)2 


a a — 1 
+ C 7 -7T2 

y [y - 1 ) 


(3.1) 


For generic choice of the coefficients a, b, c, and d solutions y = y(a) are transcen¬ 
dental functions called the sixth Painleve transcendents. These functions cannot 
be expressed by a finite number of operations through the known transcendental 
(hypergeometric-type, elliptic, etc), algebraic, and elementary functions (for a more 
precise formulation, see 

At the same time, it is known (cf. |4^) that for some particular choices of 
coefficients there exist special, and even general, solutions which can be expressed 
in terms of Gauss hypergeometric, algebraic, elementary functions (and composition 
of elliptic and hypergeometric functions). In Sect. l3.3l we show that the EFP can be 
expressed in terms of a rational solution of P6. Classification of rational solutions of 
P6 via its representation as the isomonodromy deformation of a linear 2x2 matrix 
Fuchsian ODE of the first order is given in [50]. We do not use this classification 
result, however, we need the above mentioned isomonodromy representation of P6. 

Consider the 2x2 matrix linear Fuchsian ODE 


Ay(A) = A(A)y(A), yl(A) = 4^ + ^ + ^, (3.2) 

dA AA—lA — a 

where matrices Aq, Ai, and are independent of A. In generic situation the 
isomonodromy deformations coincide with the Schlesinger deformations, and they 
are governed by the following equation with respect to a: 

^y(A) = i?(A)y(A), = (3.3) 

da A — a 

The compatibility condition of (13.2p and (13.3p . which is called the Schlesinger sys¬ 
tem, is equivalent, according to Jimbo and Miwa [4], to P6 equation (13.ip . More 
precisely, choose normalization of (13.21) such that the matrices Aq, Ai, and Aa are 
traceless, and the matrix A^o '■= — (Aq -I- Ai -|- A^) is diagonal 


A 


00 


{0oo/2 0 \ 

0 -e^/2j ’ 


doD 0. 


Denote the eigenvalues of the matrices A^ as ±0^/2, fc = 0,1, a, 00 . Actually, while 
matrices Aj- are functions of a, their eigenvalues 9k are in fact independent of a, 
since they are the first integrals of the Schlesinger system. 
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Denoting by (Afc)y, i,j = 1,2, the matrix elements of and noting the 
relation 

(^ 0)12 + (^ 1)12 + {Aa)l2 = 0, 
one can verify that the equation 


(^ 0)12 (^ 1)12 (^a)l2 

y 2/ - 1 2/ - a 


has in the generic situation (that is, for {Ai)ij + a{Aa)ij ^ 0) unique solution y, 
provided that the coefficients in m are chosen as follows: 


a = 


(0OO - 1)' 


-4 


2 ’ 


d = 


1-4 


One can associate with isomodromy deformations of (13.21) the so-called r-function, 
which is defined modulo a constant factor by the relation 


d , I'Aq + 9o/ 2 A\~\-Q\l2 -I-0^/2 

— logr = res tr ---^— -h — - 

da A=a V A A—1 A — a 


(3.4) 


= tr 


a 


A, 


a — 1 


da / 9o 01 
414"^ a- 1 


We note that our defining relation for the r-function (13.41) (with the right-hand side 
being the Hamiltonian of P6) exactly reproduce that of Jimbo and Miwa [3], since 
we use the traceless normalization of the matrices fc = 0,1, a. 

The function 


cr = a(a — 1)— log r 
da 


+ + V1V3, -I- l'2l'3)ct — 


V1V2 + V1V3 V2V3 + VlVi + V2Vi + 


where 


vi = 


0 a + 0oa 


0a — 0OO 00 + 01 ^0 — ^1 /o 

^2 =-;;-, 1^3 =-;;—, *^4 =-r—, (3.5) 


2 ’ 2 ’ " 2 
satisfies the equation (the so-called cr-form of P6): 

a^(a — l)^cr'(cr")^ -I- {(1 — 2a)(cr')^ -I- 2aa' + i^iV 2 i' 3 i^ 4 } 


= (4 + ^1) (4 + ^2) (4 + ^3) (4 + ^4) • (3-6) 


The prime denotes differentiation with respect to a. We note the cr-function can 
also be written in the form 


cr = a(a — 1)— tr 
da 


A„l + 


Ai 


a a — 1 

which will be convenient for our purposes below. 


-I- ^ir'2a — 


121122 + 123124 


(3.7) 


3.2. Riemann-Hilbert problem. Our result stated in ProDOsition l2.5l vields 
a special vector solution of the system of ODEs p.2l) . (13.31) in terms of the solution 
of the integral equation (1^ . where the matrices Aq, Ai, and Aa are given by 
(I2.2ip - (I2.23L (12.lip , and (12.121) . Here we construct a matrix fundamental solution 
corresponding to the vector one, calculate its monodromy data and, finally, refor¬ 
mulate the problem of solving of the Fredholm integral equation as the solution 
of a proper Riemann-Hilbert problem or as a solution of the inverse monodromy 
problem for ODE (13.21) . As a byproduct of our considerations here we give another 
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proof of ProDOsition l2.5l and obtain a formula for X'(A) which will appear important 
in what follows. It should be stressed that this formula, just like the statement of 
ProDOsition l2.51 can be proved by the direct differentiation, while their derivations 
based on the analytic properties of X(A) are easier modulo the considerations given 
below. 

First, we discuss the structure of the matrix V (A) in our case. Looking at (12.2111 
and taking into account that the matrix S(A) admits factorization 



we observe that the matrices Aq, Ai, and A^ are in fact residues of the logarithmic 
derivative [5 aP(A)]F(A)“^ of the matrix 

r(A) = V(A)(j W)j(l/e_(A) ,3,, 

where the functions E(A) and e_(A) are given by (12.81) and (12.6L respectively. 
Below we show that (F(A) is a fundamental solution of the Fuchs pair (13.21) . (13.311 . 
Note that the second column of y(A) is nothing but our vector solution /(A). 

The standard approach of proving that T(A) is the solution of the Fuchs 
pair (13.211 . (13.31) consists in inspecting its analytic properties. Equation (13.91) defines 
the matrix T(A) as analytic function on C \ {Cq U Coo}, i-e., piecewise analytic in 
C, with singular points at A = 0,a, 1, and oo. In case some numbers r, s, or q are 
odd one has to make one or two cuts in C connecting the singular points. If so, 
then one should consider C with the cuts, that makes no substantial changes in the 
following. 

The boundary values of the functions X(A) and E{A) suffer jumps on the 
oriented contours Cq and Coo (recall that both contours are counter-clockwise ori¬ 
ented with respect to the origin). These jumps can be found with the help of 
the Sokhotski-Plemelj formulas. Namely, let F be an oriented contour; denoting 
Yji*^(A) = F(A ± e) as e —>■ 0 with A G F, such that the point A -I- e (respectively, 
A — e) lies in the domain on the left-hand (right-hand) side of the oriented contour 
F, we have 

yr+(A) =Pr”(A)Gr, A G F, r = {Co,Coo}, (3.10) 

where 

Gco = (_\ J), = !)■ (3-11) 

Using equations (13.91) . (j2.11l) . (I2.12|) . and (12.81) one confirms the following asymp¬ 
totic behavior at the singular points: 


r(A) = + 

K k^l 





/2 


0 ^ 

(A-,z)-C/2^ ’ 


!/ = 0,a, 1, (3.12) 

where det $ 0 (^ 2 ) = 1 and the normalization condition at the point A = 00 is 

nA),;„|A + |:taoo)«„ = -(,+,). (3.13) 

The series in (j3.12l) and (j3.13l) are convergent in some neighborhoods of the corre¬ 
sponding singular points, and 9Cs are given by (12.221) . 
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We can summarize the analytic properties of Y (A) as the following singular 
Riemann-Hilbert problem: 

(1) The 2x2 matrix function F(A) is analytic in C\ {CoUCoo}, where circles 
Co and Coo are centered at 0 and oo, respectively; 

(2) The boundary values on the contours Co and Coo satisfy the jump condi¬ 
tions dSini), (IXTB : 

(3) The function Y{X) has singular points at A = 0, a, l,oo where it has 
asymptotic expansions p.l2|) and (13.131) : 

(4) The points A = a and A = 1 belong to the annulus domain between the 


circles Cq and Coo- 


The existence of the solution of the Riemann-Hilbert problem can be established 
by the reference to the integral equation (12.91) and representation (13.91) for Y (A); the 
uniqueness can be proven by standard arguments based on the Liouville theorem. 
In fact, the Riemann-Hilbert problem is a reformulation of the analytic properties 
of the matrix-function Y (A), which thus can be seen as its implicit solution. Note 
that this Riemann-Hilbert problem can solved explicitly with the help of the theory 
Schlesinger transformations. 

The solution of the Riemann-Hilbert problem delivers the generalized solution 
of the Fuchs pair (IX^ . (E31)- This can be established in the standard way by 
noticing that the logarithmic derivatives 9 aF"(A) F(A)“^ and 9 qF(A)F(A)”^ are 
rational functions on C with the first order poles at 0, a, 1, and oo. Calculation 
of the residues at these points with the help of representation (13.91) provides us an 
alternative proof of Proposition 12.51 since our vector solution /(A) is nothing but 
the second column of F(A). Note that, if one does not assume representation (13.9p . 
then the residue-matrices z/ = 0, a, 1, can be calculated in terms of the entries 

the matrix $ 1 ( 00 ). 

Using the solution of the Riemann-Hilbert problem we can construct the clas¬ 
sical solution to the pair (13.21) . (13.31) just by making analytic continuation on the 
whole complex plane of any piece of the function Y{X); since we have three pieces 
we obtain three fundamental solutions related to each other by constant matrices. 
We take the canonical solution, i.e., the one obtained by the analytic continuation 
from the disk centered at infinity with the asymptotics (13.131) . Our purpose is to 
characterize this solution by its monodromy data. If we denote this solution as 
Too (A) then it can be characterized by the following asymptotic expansions. 

Proposition 3.1. The matrix-function Yoo{X) is the unique fundamental so¬ 
lution of (13.21) with the residue matrices Ay defined by (12.211) and (12.221) with the 
following expansions at the regular singular points: 



= 0,a, 1, (3.14) 


with the normalization condition 
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where 0 , 


= —(s + q), det$o(j^) = 1 ; o.nd the matrices C^, = 0, l,a, are 

Co = GclGcl, Co^=Ci = Gcl. 

Remark 3.1. It is clear that the monodromy matrices 

My = i/ = 0,a, 1,00, 


do not uniquely define yoo(A). So, in this specific case we call the monodromy data 
Cy. Note that because integers Oq > 0, 0i > 0, and 0a < 0, the matrices Cy are 
defined by asymptotic expansions (13.141) up to the left multiplication: 


Cq —>■ 


f ao 

Vo 


bo \ 

I/aoJ 





(3.15) 


where Oy, &y G C, ;/ = 0,1, a. 


Note that our singular matrix Riemann-Hilbert problem on two contours can 
be reformulated in various ways as different singular and regular Riemann-Hilbert 
problems. In particular, one can formulate the Riemann-Hilbert problem on an 
oriented simple closed contour 7 on C, which surrounds a domain containing the 
points A = 0 and A = a but excluding the point A = 1, and with the jump matrix 
given by G-y = GcqGc,^- This can be done in two ways: first, by performing ana¬ 
lytic continuation of the solutions Too (A) and ho (A) (the solution which is analytic 
inside the disc surrounded by the contour Cq), and, second, just by observing that 
in the Fredholm kernel K one can merge the contours Co and Coo into one contour 
7 and then repeat the whole construction of this section literally. The equivalence 
of these two procedures can be established by using (13.151) . 

As a byproduct of our discussion above, one immediately obtains, with the help 
of (13.9L the following formula 

Tx(A) = .4(A)A'(A) - A'(A) (^5(0) + ■ <*^16) 

This formula is crucial in derivation of relation (13.191) below. Note that (13.161) also 
be derived directly, just by using the expressions for the derivatives of functions 
e±(A) and /±(A) entering the definition of the matrix Ai(A). Then expression (13.161) 
yields an alternative prove of the system (USD, (ESI, by noticing that the second 
term would not appear if instead of Ai(A) there would stand the matrix Ai(A) times 
some proper matrix factor, multiplied from the right. Taking into account that the 
matrix S'(A) given by (12.231) admits factorization p. 8 l) . it is easy to see that such a 
matrix factor can be chosen as in the expression (13.9p for the matrix Y (A). 


3.3. The Fredholm determinant as the cr-function. To identify connec¬ 
tion of the Fredholm determinant (12.41) with P 6 , we consider the logarithmic deriv¬ 
ative of the determinant with respect to the parameter a, 


4— logDet il — K \ = — Tr 
da V / 


( 


1 - K 


-1 



= - Tr 




(l + AjT/R.rt 


{p) Ap. 
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To evaluate these integrals, we first note that the a-derivative of the kernel K{X, /r), 
see (I2.27L can be written in terms of the functions g±{X;a), defined in (j2.17l) . as 
follows 


—K{X,g) = -^[g+{X-,a)g-{n-,a) 


47ri 


+ g-{X]a)g+{fj.]a) - 2E{a)g-{X;a)g-{fi;a)]. (3.17) 

Introducing a notation (see also (12.171) 1 


g*{X]v) = 


-5+(A;i^) 


we note that the linear combination of the functions standing in the brackets in 
(I3.17P can also be written as 

g) = -^gl (m; a)S{a)g{X-, a) 

da 47ri 

= -^ tr [5(A; a)gl(g,] a)5'(a)] , 

where the matrix S{X) is given in (12.231) . Then, using the result of action of the 
operator 1 + iZ on the functions 5 ±(A; v) in its vector form (I2.18p . we have 

d 


(n.R)_A-(.,rt 


(A) = -- 


tr 


47ri(/i — a) 

Setting here X = g and integrating over g gives 


a: {a)f{X;a)g^{g]a)S{a) 


i_ {f+{^J■)e-{g) -f+{g)e+{g)\ ^ 


dA li?—(A) -i7+-(A) 


Co M- ^^’°'^27ri fco (m “ -f-ig)e+{g)J 2ni 

d fH_+{X) -H++{X)\ 

X—a 




X—a 


Hence, we obtain the following expression for logarithmic derivative of the Fredholm 
determinant: 


log Det ("l — Tf) 
da V / 


r 

~ ~2 


^(a)X-i(a)-X(A) 


(3.18) 


A=o 


Using the A-derivative equation for the formal matrix solution, one can express 
the matrix Ar'(A) in terms of matrices X(A) and A(A), and evaluate the limit A —>■ a 
for the trace. In doing so we use the fact that even though the matrix Ar'(A) has a 
pole at A = a, the trace in (13.181) does not contain the (A — a)~^ term. To see this, 
it is useful to rewrite ()3.16p in somewhat more explicit form as 






z/=0,l,a 


2{X-v) 


{X[v)S{v)X-^{v)X{X) - X{X)S{v)) . 


Keeping A arbitrary, one can see that the v = a term does not contribute to the 
trace, since 

tr [S'(a)X-i(a) {X{a)S{a)X-^{a)X{X) - X(A)S'(a))] 

= tr [S‘^{a)X-\a)X{X) - 5'(a)X-i(a)X(A)5'(a)] = 0. 
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Setting A = a for the remaining two terms and writing 0 q,/ 2 for the pre-factor in 
(I3.18p . we have 


da 


log Det (1 — 


(i - li) 


''LX , 


= tr 


>„-i, JA-'Tia) «oX(a)5(0) , A,X(a) l),X(a)S(l)-\ 
S(a)X (a) I— -+ - 2(»-l) j 




Ai 

a — 1 


a 


-LX , 


S(a) 


0oS(O) ^ 0iS(l) 
a a — 1 


Taking into account that the matrix S(A) is upper triangular, see (12.231) . we obtain 


da 


log Det (1 — 


(i-k) 


= tr 


A, 


An 


Ai 

a — 1 


Ory (01 


01 


a — 1 


(3.19) 


Using (13.Sp . thus we have for the trace 


tr 




Ai 

a — 1 


A log Det (^1 - 7^^ + 


{i/i + t^2)(-2t^3Q + V3 + i'a) 
2 a(a — 1) ’ 


where t'fe’s are dehned in (13.51) . 

Substituting the expression for the trace above into (IrH) . we arrive at the 
following expression for the cr-function 


(T = a(a — 1)^ log Det ^1 — + (yiV 2 — viv^ — V 2 ^ 3 ) 


a 


V1V2 + V3VA - {Vl + V2){V3 + V4) 


Taking into account that from our previous discussion it follows that 

00 = s, 01 = r + g, 0 a = —f, 0 OC = —{00 + 01 + 0 a) = —{s + q), 


the life’s are therefore identified to be 
r + q + s 


Vl=V3 = — 


V2= - 


r — q — s 


Vi = 


r + q — s 


(3.20) 


Using these values, for the cr-function we obtain the following expression 

/ IN d , ^ *\ (r-l-g-l-s)^ (r-I-g-I-s)g-I-2rs 

cr = a(a - 1)—log Det [1- K] - - -- —a + - - - --. 

da V / 4 4 

It satisfies (13.61) with the v^s as specified in (13.201) . Since the EFP is represented 

in as the Fredholm determinant, Fr^s,q = Det(l — K), we conclude that we arrive 

at the statement of Theorem ll.il which is thus proved. 


4. Asymptotic expansions of the EFP in the thermodynamic limit 

Under the thermodynamic limit of a physical quantity given on a lattice one 
usually understands behavior of this quantity as the size of the lattice becomes large. 
Specihcally, in the case of the EFP besides the size of the lattice, N = r+s+q, there 
are geometric parameters, s and g, which describe the frozen rectangle at the top 
left corner of the lattice. Therefore it is interesting to consider the thermodynamic 
limit A —^ 00 preserving the geometry of the problem, i.e., keeping the ratios 
s/N and q/N fixed. Here we consider solution of this problem in the important 
particular case where g = 0; it will be convenient for us to fix the ratio s/r and 
perform the limit s —>■ 00 . We show that this problem can be solved in a systematic 
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way using the cr-form of P 6 and the asymptotic behavior of Fr^s,q at the singular 
points of this equation. 

We begin with discussion of physical interpretation of leading terms of as¬ 
ymptotic expansions in two different regimes. The details of calculations of the 
asymptotic behavior of Ti-.s.ij in each regime follow next. In Appendices we give an 
alternative derivation of the asymptotics in one of the regimes. 

4.1. Two asymptotic regimes. Recall that the EFP is defined as a prob¬ 
ability of observing the certain arrow configuration in the six-vertex model with 
DWBC. It can be seen as a ratio of two partition functions: TV,s,(j = 
where Zr^s,q and Zj\j are the partition functions of the model on the lattice with 
the frozen corner, and on the original, unmodified lattice, respectively. Since both 
partition functions in the thermodynamic limit behave as exp(—P/) where V is the 
volume (the number of lattice sites) and / is the free-energy per site, we conclude 
that in the limit, where r, s, g —>■ oo, we should have 

= exp -b o(s^)} . (4.1) 

The quantity ip is the function of the ratios r/s and q/s, and the parameter a; 
moreover, tp > 0 since Z^-gq < Zn- This function has the meaning of the change 
of the free-energy per site due to the freezing of a macroscopically large rectangle 
at the corner of the lattice. It is clear that if s > r, i.e., the bottom-right corner of 
the frozen rectangle lies below the counter-diagonal of the square, then, as follows 
from the ice-rule, Fr.^s,q = 0 and hence (p = oo. At s = r, it can be shown that for 
generic weights Fr^s,q = Zj-Zj-^q/Zj^; in our case of the free-fermion 

weights (HU = 1 for all n S N, that follows from the results of papers [HHH], 
hence Ar.s.g = (1 — At s = 0, obviously, Ar.s,g = 1, hence ip = 0. Thus the 

problem consists in studying the EFP in the thermodynamic limit for 0 < s < r. 

It is clear that ip is monotonically non-decreasing function for s S (0,r). Thus, 
there may exist a region of values of s, adjacent to zero, where ip = 0. In this case 
the quantity o(s^) in (14.11) decays much faster than s^, namely, it is exponentially 
small as s —>■ oo. In this region the EFP is close to 1, and we have the following 
behavior 

■F'r.s.q = 1 - exp{-xs-b o(s)} , X > 0. (4.2) 

Note that this kind of behavior can be expected if we recall that the EFP can 
be treated as an s-point quantum correlation function |44) and that ferroelectric 
correlations in the six-vertex model are exponentially decaying. 

It is to be also mentioned that our EFP has the meaning of a cumulative 
distribution function in the context of the random growth model, considered by 
Johansson in [34) . who proved the existence of the two asymptotic regimes dSD 
and (SU). The functions ip and x are called there as the lower and upper tail 
deviation functions, respectively. 

The existence of the two asymptotic regimes gu and (14.2p is a manifestation 
of the spacial separation of phases of the six-vertex model with DWBC in the 
thermodynamic limit, the so-called arctic ellipse phenomenon. Exactly the same 
phenomenon is often mentioned in the literature in relation to the dimer problem 
on an Aztec diamond graph (or domino tilings of an Aztec diamond) |281l5Ij . 

The phenomenon consists in a peculiar form of statistically dominating config¬ 
urations of the six-vertex model with DWBC on the N x N lattice as —>■ oo. The 
arctic ellipse arises when N x N lattice is scaled to a square of unit side length. 
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(a) (b) 

Figure 3. The Nx N lattice scaled to the unit square in the large 
N limit, the frozen square of side length ^ ^ and the 

arctic ellipse = 1 separating the phases: (a) For 

V € (it, 1) the frozen square overlaps the interior of the ellipse, in 
this case the EFP is in the disordered regime; (b) For v € (0,m) 
the frozen square lies completely outside of the ellipse, in this case 
the EFP is in the ordered regime. 


and it divides this square on two types of domains: the interior of the ellipse and 
the four corner domains, see Fig. |31 The local states (arrows) are disordered in the 
interior of the ellipse, while they are ordered in the each of the four corner domains 
due to strong influence of the boundary conditions. Recall that the EFP describes 
the probability that all arrows are ordered (frozen) in the rectangle at the top left 
corner of the lattice. There are two principal cases, whether the rectangle intersects 
the arctic ellipse or not. When the rectangle intersects the arctic ellipse one may 
expect that the EFP behaves as given by dSD. Since in this case the rectangle 
contains a portion of the domain of the disorder, we will say that then the EFP 
is in the disordered regime. If the rectangle does not intersects the arctic ellipse 
then one may expect that the EFP behaves as given by 621). Since in this case 
the rectangle encloses only a region of order we will say that the EFP is in the 
ordered regime. Indeed, when the rectangle does not intersects the arctic ellipse 
it does not essentially influence on the number of the states of the system; thus 
the EFP is close to I for large N, and, in fact, equals to I in the thermodynamic 
limit. Contrary, when it intersects the arctic ellipse it starts to order the disordered 
domain and reduces significantly the number of states, so that the EFP should be 
small for large N, and, in fact, vanishes in the limit. In Fig.[3]these two situations 
are shown in the case q = 0, i.e., when the rectangle is just square. 

In calculations we limit ourselves here to the case q = 0. In this case the 
analytic results appear considerably simpler without altering the main ingredients 
of the method. We take the two remaining integer parameters, s and r, 0 < s < r 
to be large in the limit, with their ratio being a non-vanishing parameter 

v=-, i;e(0,1). 

r 


24 
















(4.3) 


The condition of intersection of the square with the arctic ellipse reads 

1 - \foi 

V = U, U= - 

1 + Va 

For V G {u, 1) the EFP is in the disordered regime, and for ?;(0, u) it is in the ordered 
one. Sometimes instead of natural variables v and u, it is useful to use the variables 
a and /3, in which the condition of intersection of the ellipse reads 

a = /3, P= 

For a € (/3,1) the EFP is in the disordered regime, and for a G (0,/3) it is in the 
ordered one. 


4.2. Asymptotic expansion in the disordered regime. We begin with 
the study of the EFP in the disordered regime, which corresponds to n G (u, 1) or 
a G (/3,1), by considering behavior of Fr^s,q at the critical point a = 1 of P 6 . 

It was noticed in [3^ that C'r,s(l — where the constant Cr,s is 

essentially given by the value of the Hankel determinant in (12.21) at a = 1. However, 
this asymptotics does not uniquely characterize the corresponding solution of the 
(T-form of P 6 , and therefore we need a more elaborated result. 

Proposition 4.1. Ts a —>■ 1, the EFP behaves as follows 


Fr,s,o = a.s(l - ay‘ ( 1 - - a) 


s(r — s) (2s^r — 2s^ — 3s^ + l) , „ \ ^ ^ 

+ - 4(tf-l) -i(l-a)VO{(l-a)>)j. (4,5) 

where the constant Cr,s has the explicit form 


s-l 


Cr,s = n 




j=0 




(4.6) 


Proof. At a = 1 the Hankel determinant in (12.21) appears to be related in the 
canonical way with an ensemble of Hahn polynomials (SOllHl]- At g = 0 we deal 
with the special case of the Hahn polynomials Qn{x;0,0,r — 1); for the notation 
and basic properties see, e.g., [53]. We will use the normalized version of these 
polynomials with the coefficient of the highest-order term equal to 1 , 

’ —n, n + 1, —X 


Pn(^) — / I i\ 1) — I 1 \ 3-^2 


1 , -r + 1 


(R + l)n " ^ ^ (n + 1 ). 

where (a)„ = a{a+ 1) • • • (a + n— 1) is the Pochhammer symbol. These polynomials 
satisfy the orthogonality condition 

r — 1 


^ '^ Pn{F)Pn'{F) — hjiSnn'^ 


hn — 


(n!)^(r + n)! 


x—Q 

and the recurrence relation 
r —1 


( 2 n)!( 2 n + l)!(r — n — 1 )! 

n^(r^ — nf) 


hn 


(4.7) 


■ (4.8) 


XPn{x) = Pn+l{x) + ^—Pn{x) + - - Pn-l{x), ^ A(A 2 

2 hn-1 hn-1 4(4n"^ - 1) 

Let us denote the Hankel matrix in the representation (12.21) by FI (a); we also write 
for brevity H = H{1). At g = 0 (for the case of g 7 ^ 0, see |S5]) the entries of 
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the matrix H appear to be the moments of the orthogonality measure in (HIT) and 
therefore 


detH = det 

l<j,k<s 


jn—{) 


s—1 s — 1 

=n =n 

j=0 j=0 


(2j)!(2j + l)!(r-j-l)!- 


Taking into account the value of the prefactor in (12.21) . we obtain the expression 
(ITO for the constant Cr^s- 

Let us now focus on the terms of the (1 — Q;)-expansion in (14.51) . We de¬ 
note H' = ^H{a)\a=i and H” = ■^H{a)\a=i. Using the fact that detiL(a) = 
exp(tr log iL(a)) and differentiating, one can show that the coefficients of the Taylor 
expansion 

det H{a) = det H ■ -|- ci(l — a) -I- ^(1 — -|- O ((1 — a)^)^ , 

can be written as 


Cl = - tr(iJ'iJ-i), C2 = {tr{H'H-^)f - tr{H' -h tr{H"H-^). 

Now we use the fact that 

„-i ^ 1 1 ^ pn{x)pn{y) 

Taking into account that “ !)> 

one can reduce evaluation of the traces in the expressions above to the standard 
sums involving the related polynomials, namely 

mpl{m) _ r -1 ^ w?pi{m) _ hn+i , 

, hn 2 ’ ^ hn hn \ 2 J hn-1 ’ 

which can be easily derived from the orthogonality condition Kn and the recur¬ 
rence relation (14.8|) . In this way, for the coefficient Ci we get 


Cl 


= -EE 


mp^{m) s(r — 1) 


n n 2 

n—0 m—0 

For the coefficient C 2 , after the straightforward but more involved calculations, we 
obtain 

s(r — 1) s^(l — 2r-I-s^) — s^) 


C2 = — 


4s2 - 1 


Finally, taking into account that 




q,s(s-1)/2 ■ 2 

we arrive at (|i3]) . 


□ 


As a direct consequence of (14.51) . we have the following. 


Corollary 4.1. As a —>• 1, 


logFr,s,o = logCr.s -b S^ log(l - o) - ^E-E(l - a) 


s{r — s) (7s^ — sr — 2) 


8(4s2 - 1) 
26 


(l-a)2 + 0((l-a)3) . (4.9) 
























To construct the asymptotics of logFr^s,o for large s and r = s/v, we need the 
corresponding behavior of the coefficients of (14.91) . A nontrivial calculation concerns 
only log C'r,s- To do that we rewrite (14.61) with the help of the Barnes G-function 
[54l as follows: 


Cr.s = 


^s+ 1/2 G'^{l/2)G{r + s + l)G(r - s + 1) 
2A2 ^-i) G2(r + l)G(s + l/2)G(s + 3/2) ‘ 

The Barnes G-function satisfy the relations 

G{z + l) = T{z)G{z), G(l) = G(2) = G(3). 

As established by Barnes: 


(4.10) 


logG( 2 ; + 1) = y logx- jz^ + -;^logz + C'(-l) 


log27r 


1 


12 




B2k+2 


2k 


^ 4fc(fc + 1) 


+ (4.11) 


where = —0.165142... is the derivative of the Riemann function ((z) at 

z = — 1, and i? 2 fc’s are the Bernoulli numbers: 


Bo-1, i? 2 -y 


Bo = 


1 

42’ 


B. = -, 


1 


Bio = 


66 ’ 


(4.12) 

Therefore, taking the logarithm of (I4.10p and using (14.111) . we arrive at the 
following result. 

Proposition 4.2. As s,r —>• oo, with s/r = v fixed, v G [0,1], we have 
log Gr,s = -s^ {log 4u - log(l + u) - log(l - v) 


1 


2v^ 

1 (1 — u^) 


-log s-log 

12 ® 12 ® 


where the coefficients B 2 k = B 2 k{v) are 


2v^ 

+ c'(-i) + E% + o(^ 

k=i ^ 


% (4.13) 


B2k = 


,,2k 


.,2k 


B2k+2 

4fc(fc + 1) \^(l+u)^^ ' (1 — u)^^ 

fc-1 

-E 


- 2v 


2k 


R. 


2{k—m)-\-2 


m—Q 


2k- 1 


^ 22”i+i(fc — rn)(fc — m + 1) \ 2m J 2'^^+^3k{k + l){2k + 1)^ 


4fc^ + 6fc — 1 


where (^ 2 m^) Binomial coefficients. For example, 


B2 = —. 


1/1 1 + 


15 


^^(l + u^) 
15(1 -U2)2 


, Bi = - 


u®(u®-4v‘‘ + 5w2-10) 31 


504(1-z;2)4 


16128 


Note that, besides single logarithmic term, the expansion for log Gr,s involves 
only even powers of 1/s. The same property appears to be valid for the explicit 
terms in gH). 
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To derive expansions of the EFP in the thermodynamic limit, we start with 
analyzing the cr-form of P6 in the large s limit. We note that for q = 0 the a- 
function in (EH) reads 

Absence of the odd powers of s in the cr-form of P6 (ll.3p together with the expan¬ 
sions (14.911 and (14.1311 suggests that the cr-function may be searched in the form of 
the following asymptotic ansatz in the decaying powers of s^: 

a = CT 2 S^ + ao + <J-is~'^ + • • • j s —t c». (4-15) 



Note that if we succeed in the construction of the expansion (14.1511 with the coeffi¬ 
cients CT 2 fc = cr 2 fc(a), which are analytical functions of a, then the Wazow theorem, 
see Th. 36.1 in Chap. IX of [55], implies that there exists a genuine solution with 
asymptotics (I4.15L To justify that this solution actually coincides with the solu¬ 
tion given by (14.1411 one can verify that both solutions have the same asymptotics 
as a —>■ 1. Here we present the general scheme of derivation of asymptotics; more 
details are given in Remark 14.31 below. 

The expansion (14.1511 can be constructed in a standard way by substituting it in 
EH- On this way we first obtain the leading term, by requiring that it reproduces 
the asymptotic condition at a —>■ 1, and next prove that the further terms can be 
obtained recursively. 

It happens that the leading term a 2 can be found as a solution of the first- 
order ODE which is obtained by noticing that the first term in the left-hand side 
of (11.311 . which contains second-order derivatives, is sub-leading as s large. Besides 
the overall factor tr^, the resulting first-order ODE splits on two equations 


1 


11 — /1 -(- 1)“^ \ 

+ CT2 -b (1 - a)cr2 = 0, —-^ ( ^^^2-^2 + aa'A = 0. (4.16) 


4^2 ■ ■ V - -;-2 

These equations have the following general solutions 


(o’ 2 )i = Ci (a — 1) — 


4w2 


(o'2)ii = Cii a -b 


1-b(1 — 


2\2 


4v^ WCiiv'^ 


(4.17) 


respectively. Here Ci and Cn are integration constants, which may depend on v. 
Moreover, the second equation has two partial solutions 

1 -b j;2 1 — 7;2 

('^2)± = ± (4.18) 

Consider now the function (p, defined by (0311. Using (14.1411 and (I4.15|l we find 


that 


ip = 



(1 -b v)'^ 
4^2 


a — 



da 

a(l — a) 


+ C, 


where C = C{v) is an integration constant. The function tp should obey the fol¬ 
lowing asymptotic condition as a —> 1: 


tp = 


(1 -b v']^ (1 — rj)^ 

- log(l - a) -b log4?;-log(l + v) -“ '^) 


1 — V 
2v 


(1 - a) - 


2^2 

(1 - u)(l - 7v) 


32z;2 


{l-af+ 0{{l-af) . (4.19) 
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The function which corresponds to solution (< 72)1 reads 

This function does not have asymptotics (I4.19|) , because of the absence the leading 
term — log(l — a). 

In the case of the function (<p)ii, which corresponds to (o' 2 )ii, the leading term 
— log(l — a) of the expansion (j4.19|l dictates the choice Cu = — then we get 

1 — V 

(<p)ii = - log(l -a) --— log a + C. 

2v 

As (X —y 1, 

((/3)ii = - log(l -a) + C+ ^-^(1 - a) + -^^(1 -af + 0 ((1 - af) . 

Therefore, because of the coefficient of the term (1 — a)^, the asymptotics (14.191) 
cannot be satisfied for any choice of C. 

Now consider the functions corresponding to partial solutions (I4.18|) . 

(‘^)± = - log (1 ± \/a) - log (1 T Va) + log a + C. 

Asymptotics a —>■ 1 of the function ((/3)+ has the wrong coefficient of the leading 
term, namely, — 1/u^, compare with (14.191) . For the function {(p)-, as a —>■ 1, we 
get 

{p)- = - log(l -a) + C - ^ log 2 

Therefore, comparing the last asymptotics with (14.191) . we observe that (p is given 
by the function (p)- with the proper choice of the integration constant C, 


V = - log j - ^ log j + log a 

(1 + (1 — 

+ logu-^-^log(l + u)-^-^log(l-u). (4.20) 

Note that p is positive for a G (/3,1), because it can be presented in the form 



where /3 is given by (lOl) . 

Formula (14.201) . when written in terms of the parameter u, is exactly the ex- 
pression (fra and coincides with the leading term of the asymptotic expansion of 
the EFP in the disordered regime obtained in m- 

Thus we have just constructed the term 172 of the expansion (j4.15|) ; the further 
terms a- 2 k, fc = 0,1,..., satisfy the recurrence relation of the following form 

[v^{l + - (1 - a-2k = P{v,a;a2, ■ - . ,cr_2fc+2), (4.22) 
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where P is a polynomial in all its variables. Thus all coefficients (t_ 2 /c can be 
recursively constructed as analytical functions of a G C defined in the disc centered 
at a = 1 and of radius 1 — /3. 


Remark 4.1. Similarly, one can construct expansions of the form (14.151) which 
correspond to remaining functions in (14.171) and (14.181) . namely, to (<^ 2 ) 1 , (<^ 2 ) 11 ; 
and {jp 2 )+- According to the Wazow theorem they correspond to different (more 
general) solutions of the a-form of P6. 


Remark 4.2. The general solution of the a-form of P6 JOl) with q = 0 is 
determined by two parameters in the expansion of the form (Ol) . namely, by the 
coefficients of the terms log(l — a) and (1 — a)^. Our construction shows that if the 
coefficient o/log(l —a) is fixed as then all solutions whose coefficients of — 
as s —>■ 00 and s/r = const., are of the form (r —s)(7s —r)/32 + o(s^) have the same 
asymptotic expansion (|4.15|) with the leading term given by (o' 2 )-. Therefore their 
asymptotics as s —>■ 00 differ by exponentially small corrections which are defined 
by the function of o{s^). To find these corrections seems an interesting problem. 


From (|4.14p and (14.151) it follows that 


log Fr,s ,0 = -s^</5 - — log s 


E 


Q2fc 

c2k 


0 (s 


where the coefficients a 2 k = 02 ^ (a, u) can be obtained by integration of the cor¬ 
responding coefficients a- 2 k and the integration constants are fixed by expansion 
(I4.13p . For ao we find 


1 / 4v^^/a \ 1 1 — 

^ 8 (,2(l + u2)y5-(l-u2)(l + a) j " 12 

Similarly, one can obtain further terms; for example, we find 


C'(-I)- 


02 


and 


(1 - v^)(l - ^/a)^ [ 2(1 -I- -I- 5(1 - a) 


2 2 

V 


64 [2(1 -h - (1 - v2)(l -h a)]" 


1 /I 


(1 - 


1 -I- v^(l -I- v^) \ 

15 15(1-u2)2 J 


(24 = 


(1 - u2)(l - 75)2 


8(1 -I- 7a)®(l — 375 -I- a)' 


,10 


256 [2(1+ v^)^- (l-i;2)(l + a)] 

-I- 10 (1 — 7 ^^ (1 + (5 — 467a -h 5a) 

- 5(1 - a)^ (11 -b 10675 -b 11a) u® 


-b (1 — 75 )^ (1 -b 75 )^ (1 -b 187a -b a) v'^ 

— (1 — 75 )* (5 -b 1475 -b 5a) -b (1 — 7a)^° 

u®(u® - 47 -b 57 - 10) ^ 31 

504(1-7)4 16128’ 


Calculations show that the functions a 2 k become more and more cumbersome as k 
increases. Finally, rewriting the above formulas in terms of the parameter u instead 
of a, see (H31) . we arrive at Theorem o 
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Remark 4.3. Our justification of ansatz (I4.15P in fact includes the following 
steps: 

(1) Appearance of the formal expansion (14.151) as the simplest possible analytic 
(in a) generalization of the behavior of Fr^s,o for large r, s (s/r = const. = 
v) at a = 1. 

(2) Substitution of the ansatz into the a-form of P6 equation to prove that 
this ansatz solves the equation. 

(3) Application of the Wasow theorem /' |55j , Th. 36.1) to prove that there 
exists at least one genuine solution of P6 that has asymptotics (14.151) . 

(4) Proof that the solution obtained via the Wasow theorem coincides with the 
function a related to fr.s.Oj see (|4.14p . based on matching of asymptotics 
expansions as a ^ 1 using the double character of asymptotics (I4.15p with 
respect to s ^ oo and a —>■ 1. 

Some of the items of Remark 14.31 require further comments. 

Concerning item (2), using mathematical induction and formulas (I4.18p and 
(I4.22P one can prove that the coefficients a- 2 k, k = 0,1,..., are rational functions 
of ■\/a with the poles only at the point a = /3. Obviously, these functions are 
holomorphic with respect to a in any simply-connected domain in C \ {0, /3}. 

Concerning item (3), the Wasow theorem does not apply to the cr-form of P6 
directly, because it deals with the first-order vector ODEs resolved with respect to 
the derivatives, with holomorphic right-hand sides, and a small parameter e. This 
theorem applies instead to the corresponding representation of P6 as a Hamilton¬ 
ian system. The cr-function is intimately related to the Hamiltonian, and is given 
as a quintic polynomial in terms of the canonical variables [46]. Conversely, the 
canonical variables can be written as rational functions of cr, da<x, d\a and a (see, 
e.g.. Sect. 2.1 in [46]). Therefore, for the canonical variables of the Hamiltonian 
system there exist formal power-like expansions in e, where e = 1/s. The Wasow 
theorem also deals with two assumptions, (A) and (B). The assumption (A) con¬ 
cerns basically the definition of the domain (in the variable of) of validity of the 
asymptotics; in our case it is evident that it is nonempty, and contains any segment 
in the interval (/3,1). The assumption (B) concerns nongeneracy of the matrix of 
the vector ODE, after, possibly, a proper transformation of the canonical variables. 
In our case it follows from the fact that our e-expansions of the canonical variables 
are defined uniquely as long as the leading term is given. An example of such 
transformation for the Hrst Painleve equation is given by Wasow (se^ |^, p. 225). 
Our case is consired in Appendix |A| 

Concerning item (4), the matching of asymptotics is made by using the known 
asymptotics as a —?> 1 up to the terms of O ((1 — a)^). It is known that this 
asymptotics defines uniquely the solution of the cr-form of P6. Therefore, by com¬ 
paring the expansions as a —>■ 1, one can verify that the solutions, a priori defined 
in a different way, actually coincide. Strictly speaking, the Wazow theorem guar¬ 
antees existence of at least one solution with the asymptotics (14.151) . The other 
possible solutions may differ by terms decaying faster that any integer power of e 
(e.g., exponentially small terms). In our case, because of exact coincidence of the 
coefficients (no free parameters) of the terms (1 — a)^, fc = 0,1, 2, the exponential 


^There is a misprint after Eq. (36.40) in 


m- the assumption (B) is meant rather than (A). 
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small corrections (if any) does not influence on the matching procedure. In other 
words, the corresponding monodromy data of the solutions coincide. 


4.3. Asymptotic expansion in the ordered regime. To fix integration 
constants for the terms of the asymptotic expansion of the EFP in the ordered 
region, v € (0,u) (or, equivalently, a € (0,/3)), we use the following result about 
the EFP in the limit a —>■ 0. 


Proposition 4.3. As a —>■ 0, the EFP behaves as 


F — 1 — 

r,s^q — 


where the standard notation for the binomial coefficients have been used, = 


/c!(n —fe)! ‘ 

Proof. Consider the n = 1 term in the sum (12.311 . which is just the trace of 
the our operator K, see Proposition 12.31 


Tr(A) = 

Change variables: 


(A — aY d\ f {w — dw 


Co ~ 27ri {w — aY{w — A)^ 27ri 


A = av, w = l//r. 

Note that since 0 < a < 1 contour Co —>■ Co, while Coo —>■ — Co, thus 


Tr(A) =a’-'*+^(-l) 


(1 - vY dv 


(1 - ^^Y+'^ 




Tco (1 ~ avY~^'^^^ 27ri Jq^ (1 — Q;/i)’'/j,®+®(l — av^Y ^Tii 

Therefore as a —>■ 0 we find 

/ (1 - uY du r (1 - /i)’'+9 d/i 
fco 27ri Jc^ 27ri 


Tr(A) =a’-"+^(-l) 


(l + 0(a)). 


Now we calculate the above integrals via residues, that gives (—and 
(—respectively, and so we get the second term in the right-hand 
side of (14.231) . 

Analogous consideration shows that Tr(A^) = 0(0;^^’'“'’+^^) as a —0 and thus 
the further terms in (12.31) do not contribute to the leading order of the Taylor series 
of Fr^s,q at a = 0. □ 


Setting g = 0, we have the following result. 
Corollary 4.2. As Q, 

log(l - Cr,s.o) = (r - s -M) log a 2 log 


s — 1 


0{a) (4.24) 


and 


log Fr,s,o = - exp (r - s -I- 1) log a -I- 2 log 
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s — 1 


+ 0(a)L (4.25) 











Our next task is to get asymptotics of the logarithm of the binomial coefficient 
in (I4.24P as s —>■ oo, r = s/v. We recall the known asymptotics of logarithm of the 
Gamma-function as z —>■ oo (see, e.g., [5S]1. 


logr(2;-ba) = [z + a-^j\ogz-z + i log(27r) 


+ f2 


n—1 


n{n l)z^ 


I argzl < TT, 


where Bn{a) are the Bernoulli polynomials, 


Bn{a) = ^ 




Bka 


i—k 


Here Bk are the Bernoulli numbers, Bi = —1/2, -B 2 fe+i =0, fc = 1,2,..., and i? 2 fc 
are given by (I4.121) . 

Using the above formula for asymptotics of logarithm of the Gamma-function 
we arrive at the following result. 


Proposition 4.4. Hs s —>■ oo and r = s/v, 


2 log 


s — 1 


-2 ( log u + log(l - s - log s + log ^^3 


E 

n—1 


- (H2„ + 2n)(^)2"-i - B 


^2n 


i(2n — 


E 


,,2fe 


^ fc(l — vY^s"^^ 


-b 0 ( 3 -"™-^). (4.26) 


Gorollarv 14.21 and Proposition 14.41 imply that the simplest possible asymptotic 
ansatz for the logarithm of as s—>- 00 ,r = s/i;, reads 


log Fr,s,o = - exp < -xs - log s -I- 


0 (, 


— m — 1 


) • 


(4.27) 


fc =0 


In view of (14.141) . the asymptotic ansatz for the cr-function, which corresponds to 

(023, is 


u = I + 1- ) »“ +e-‘" +0(6-’“"), 


4., -■2.,,"- 6 

where the function lo admits an asymptotic expansion in the inverse powers of s, 

00 = ojq U 01 S ^ -)- 002 S ^ -)- * * ■ . (4.29) 


At first glance, the asymptotic ansatz (14.281) is considerably different in compar¬ 
ison with the ansatz for the cr-function in the disordered regime, see (14.151) . How¬ 
ever, it can be viewed as the same ansatz but with the function ip = 0, a G (0,/3). 
Indeed, as it is mentioned in Remark l4.21 an asymptotic expansion for the cr-function 
may contain exponentially small corrections, which are not indicated in (|4.15l) . At 
the same time, in case ip = 0, the corresponding function 172 is just given by the 
coefficient of the first term in (14.141) . From the point of view of the solutions of 
(I4.16|l , which govern a 2 , this expansion corresponds to the unique case where both 
general solutions (14.171) coincide, Ci = Cu = —(1-|-u^)/4u^. It happens that in this 
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very case a- 2 k = 0, k = 0,1,_ Comparing (14.141) and (I4.28P with (14. 2p we see 

that the function y is equal to uq, moreover it should be positive for a € (0,/3). 

Justification of asymptotic ansatz (14.281) . (j4.29l) can be obtained analogously 
as it was done for ansatz (|4.15|) . However, here we have an independent proof for 
ansatz (14.271) based on the saddle-point method applied to the Fredholm determi¬ 
nant, presented in Appendix |B] Therefore in what follows we consider an explicit 
construction of the coefficients in the series (j4.29|) . 

From the cr-form of P 6 it follows that the w-function in all orders in 1/s obeys 
the following ODE: 




= [(a — l)sa;' -I- 1] 


(1 -I- v)^a — (1 — v)^ 


(1 -I- vy 


(4.30) 


V V 

For example, keeping the leading order in this equation, we find that the function 
X = Wo satisfies equation 

(1 -I- v^a — (1 — vy 


\a-l){^y = 


(4.31) 


Since a G (0,1) and y' is real, the right-hand side of (14.311) is negative, therefore 
a G (0,/3), where /3 < 1 is given in (14.41) . Using (I4.24|) and (I4.27L we note that y 
has the following asymptotic behavior as a —>■ 0 , 


^ ^ _ qj 

X = - logo; -I- 2 logn -I-log(l - r^) ) -|- 0{a). 


Solving (14.311) for y', we obtain 


x' = -— 

va 


(1 — vy — (1 -I- vyc 


/3 — a 


(4.32) 


(4.33) 


1 —a (1 — ^ 7 ^) 0 : V 1 — a’ 

where sign minus is chosen to ensure the correct coefficient of the log a term in 
(I4.32p . Integrating (14.331) and fixing the integration constant according to the sec¬ 
ond term in (14.321) . we obtain 


X = 


4 [ - a) + yy - a\ ( 


Vi — a -I- yy — < 


Note that this expression can also be rewritten as 

r, 1 + y/P 1 ( y/P + Q: -I- ^ 

X = — 




(1 -h yy)y/a 


- 2 log 


yy — a + \/{l — a){y — a) 

(1 - yy)ya 


, (4.34) 


or, noticing that the numerators of arguments of the logarithms are perfect squares, 
as 


1 -I- yy ( a/(i - 1 - ya){yy+ya) + \/(i — ya){yy — ya) 

X = 4:;- 7 = log I - 


i-vy 


y 2 {i + yy ) y ^ 


_ ^ / \/{I - yayyy+ya) + ^ {I + yayyy - ya) \ 

\ \/2(i - yy)ya j 


(4.35) 
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These expressions imply that x\a=p = 0, and furthermore y > 0 for a G (0,/3), 
since it can be written, similarly to (14.211) . as the integral 

Using (14.341) it can be shown that y = 2J{t + 1), where J{t) is the upper tail rate 
function (for 7 = 1) obtained by Johansson, see Eqs. (2.21) and (2.25) in [^, and 
< = — 1. R.ewritinff ()4.35D in terms of v and u. see (14.31). we arrive at dill). 

After we have determined wq = y, the first derivatives of the further coefficients 
in the expansion (14.291) . namely w(.’s, can be derived uniquely from (14.301) . Equa¬ 
tions (14.141) . (14.271) . and (14.281) implies the following equation for the coefficients 

h- 


K 


log a + log(l -a) + log(-y') -Hog 1 + ^ (-UlV+i +j 

V k=0^ ) k=o ^ k=0 


W/c+l 


For example, the first three coefficients 6 ^ are 

bo = -oJi - log(-y') - logo - log(l - a), 
u K 

hi = —UJ2 H-7, 

X 

fa--“> + 7 +2(77 

The undermined constants in are fixed with the help of the expansions (14.251) 
and (14.261) . Using this algorithm, we obtain the following expressions: 


bo = log 
bi = - 

b2 = 


y2TT\/l — a [(1 — n)^ — q:(1 -|- j 

(1 -I- n"^))! — a)^ -I- 9{v + r^^)(l — a^) — 2u^(10a^ — 11a -|- 10) 
6\/l — a [(1 — u)2 — a(l -I- 

■ [(a^ -I- 8a -I- 1)(1 - a)^(n'^ -|- 1) 


(4.36) 


(1 — a)((l — n)^ — (1 + f)^a)3 
— 4(1 — a)(a^ -I- l)(z^^ -I- 'c) -f 6(1 — a)"^ti^ -f 4a(a^ -I- a -I- l)n^]. 

Clearly, 

m j 

log(l - Fr,s,o) = -xs - log s + X! “I + 0(s“'"“^), 

k=0 ® 

that finishes the proof of Theorem 11.31 

It is interesting to note that (14.301) can be solved explicitly in terms of the 
hypergeometric function, see Appendix [Cl 


5. Conclusion 

The goal of this paper is twofold: to elaborate an algorithm for construction 
of full asymptotic expansions for a physically interesting quantity for an integrable 
model, and to show that the theory of P 6 can serve as an effective tool for this 
purpose. 


35 



































The situation is not absolutely standard since we used an ODE with respect 
to a continuous variable to find asymptotic expansions with respect to discrete 
parameters. At first glance, Toda or discrete Painleve equations might be more 
suitable for this problem. Indeed these equations can be used in obtaining the 
leading terms of our asymptotic expansions 130], however, a perspective of obtaining 
higher order terms (HOTs) of the expansions with the help of the difference or 
differential-difference equations seems unclear. 

Our problem looks like a standard object for asymptotic methods based on 
the Riemann-Hilbert conjugation problem or isomonodromy deformations (see Sec¬ 
tion IQ) . These methods however deliver the leading terms of the corresponding 
asymptotic expansions, while the explicit formulas for HOTs require much more 
deliberate and technically complicated efforts, so that anyway HOTs are easier to 
derive via substitution of the asymptotic expansions into the corresponding inte- 
grable (differential, difference, etc) equation. 

In this paper we show that in our case not only HOTs but also the leading terms 
can be obtained via the substitution of asymptotic expansions into P6. Since we use 
the ODE that plays an auxiliary role with respect to discrete parameters we need 
initial conditions for determination of all terms of the asymptotic expansion rather 
than just for the leading terms, as it usually happens in the standard approach. To 
get these initial conditions one needs to know asymptotic behavior of the r-function 
at critical points of P6. In our case it is asymptotic behavior of Fr^s,o as a —>■ 1 
and a —>■ 0, Propositions 14.11 and 14.31 respectively. The results stated in these 
propositions are, in fact, equivalent to the connection formula for the P6 r-func- 
tion. We have found no corresponding result in the literature, probably because 
our solution is too special from the point of view of the general theory of P6. 

Asymptotics as a —>■ 0 is obtained in Proposition 14.31 with the the help of the 
Fredholm kernel AT, which is the principle manifest of integrability. At the same 
time the proof of asymptotic expansion as a —>■ 1 in Proposition 14.11 is obtained 
by an ad hoc observation, which may not occur for some other analogous problem. 
We used this observation because it simplifies our derivation, however we could 
employ a more deliberate scheme, we are going to outline below, where a —>■ 1 
asymptotics is not required. In principle, the connection problem for r-function of 
P6 is a different and external problem in comparison to the one we consider here, 
so that we can just refer to it. However, there is a way to avoid this problem, 
or more precisely, to change it to another one: the connection problem for the 
“corresponding” r-function of the second Painleve equation. 

The key point of this scheme is the representation of Fr^s,q as a determinant 
of the integrable Fredholm operator K. The next step is the derivation of a-form 
of P6. Our problem is to find asymptotics of Fr^s,q for large variables r, s, and q 
in the thermodynamic limit {s/r = v = const) for all values of the parameter a, 
0 < a < 1. As shown in the paper, it is easy to find asymptotics in the ordered 
regime, i.e., in the interval 0<u<u = (l — i/a)/(l -I- y/ot). This asymptotics 
requires a —>■ 0 asymptotics and the cr-form of P6 for DetAT. In Appendix iBl 
we show that in this case we actually can find asymptotics directly, without any 
reference to P6, by the help of the saddle-point method for the Tr K. 

Then we observe that in the neighborhood u v our asymptotic expansion 
is destroyed (see Theorem Ol) and the so-called “boundary” variable z can be 
introduced, which, as is easy to guess from the structure of HOTs, reads z = 
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{u — v)s^/'^. Substituting it into the tr-form of P6 instead of a one obviously finds 
that (j(a) —?> 0 - 2 ( 2 ), where 0 - 2 ( 2 ) = f y^(z)dz is a solution of the cr-form of the 
second Painleve equation (P2). It is well-known that 0 - 2 ( 2 ) = — f y‘^(z)dz, where 
y{z) is the corresponding solution of the canonical form of P2, 

y" = 2y^ + zy + C. (5.1) 

Let us explain why the appearance of P2 is an obvious thing here. As follows 
from the calculations presented in Appendix [B] there is a coalescence of two simple 
saddle-points in Tr iL dX u = v. This is equivalent to coalescence of two simple 
saddle-points in the steepest-descent method for the corresponding RH problem or 
coalescence of four turning points in the isomonodromy deformation method. This 
process is described by the solution of the model RH-problem corresponding to P2 
(see details, say, in m)- Surely, one can deduce it directly from the corresponding 
differential equations with the help of the above-mentioned substitution and obtain 
one of the limits P6 -X P2 (see, say, discussion in introduction to [S5]). How to 
choose the proper solution of equation (EH) in the framework of this approach? 
The solution of P2 should match with the asymptotics in the ordered case: this 
means that in terms of variable 2 it should be exponentially decaying: thus the 
constant C in EH vanishes. It is well-known that there exists only one-parameter 
family of solutions of P2 regular on the positive real semi-axis and exponentially 
decaying there. The value of this parameter should be determined via the matching 
procedure in the domain 2 ^ +5*^, e > 0, with the leading term of asymptotics in 
the ordered regime, v < u, 

On the negative semi-axis solutions of this family can be oscillating or possess 
infinite number of poles, depending on the value of the parameter of the family. 
There exists also the unique solution that separate these two types of asymptotic 
behavior, it can be treated as the solution with the only pole at the point at infinity, 
it is known as the Hastings-McLeod solution [^. This solution grows as -y/— 2/2 
for 2 —>■ — 00 . 

Since we already know asymptotics of Fr,s,q in the disordered regime, which is 
regular and non-oscillating, it is obvious that the only suitable solution of P2 is the 
Hastings-McLeod solution. However from the formal point of view this fact should 
be established by applying the matching procedure described above. 

Asymptotics of the Hastings-McLeod solution for 2 ^ e > 0, serves as 
the “initial data” for asymptotics of Fr,s,q in the region v > u. Obviously, to 
use the “P2-asymptotics” as the “initial data” we have to develop this asymptotic 
expansion up to all decaying orders of 2 . From our point of view, this expansion is 
an important ingredient of the asymptotic scheme based on matching of asymptotic 
expansions. In certain cases, the scheme which can serve as an effective alternative 
to other asymptotic methods. 

At this stage we just explain that the leading terms of the “P2-asymptotics” 
and the expansion constructed in Theorem 11.21 do match. For the leading term we 
have 

\ogFr,s.q^- f (z- z)y^{z)dz-{v-u)^s^/12, (5.2) 

J —00 

since y ^ -^— 2 / 2 . This shows that the “Painleve asymptotics” in the transition 
domain matches at least with the leading order of asymptotics obtained in The¬ 
orem [L2l Above we have given an explanation of appearance of P2 purely in 

37 



terms of the theory of differential equations and nonlinear special functions. For 
the corresponding matrix model the result is well-known via the appearance of the 
“Airy-kernel” of the associated integrable integral Fredholm operator K [60]. 

Simultaneously, the last asymptotic relation in (15.21) explains the physically 
interesting phenomenon of the third-order phase transition mentioned in Intro¬ 
duction: the leading term of the asymptotic expansion of log FV,s,g (denoted as (p 
in Theorem HH) vanishes at u = w, together with the first two derivatives (with 
respect to v or u), whilst the third derivative does not. In other words, inter¬ 
preting Lp as the change of the free energy per volume in the thermodynamic limit, 
—v"^p/{l + v)'^ = lims_r._).oo \-OgFr,s,q/{r + sY, one concludes that (j5.2l) implies p = Q 
for V < u, thus reproducing the result of [30] about the leading term of asymptotics 
described by the function p. 


Appendix A. The Wasow theorem and the cr-form of P6 

The cr-forms of Painleve equations does not fit the setup of the Wasow theorem 
[Th. 36.1 in Chap. IX of |5S]], because they are quadratic with respect to the 
second derivative of the cr-function. Nevertheless, in many cases this theorem can 
be applied to these equations due to the presence of the Hamiltonian structure. 

Here we give some more details to the general comments given at the end 
of Sect. 4.2 that support item (3) of Remark 4.3. Namely, we show here that 
the Wasow theorem applies to the representation of P 6 as a Hamiltonian system, 
and due to the one-to-one correspondence between the canonical variables and cr- 
function (as proven by Okamoto |46j 1. it implies the validity of the similar result 
for the cr-form of P 6 . We first consider the P 6 as a Hamiltonian system in a general 
setup, and next proceed to our case. 

Since here we deal with P 6 as a Hamiltonian system, we will rely upon results 
obtained by Okamoto [46j . To keep contact with our discussion in the main text, 
we make the following identification between parameters vi,... ,Vi of Jimbo-Miwa 
and those 6 i,..., &4 of Okamoto: 

— ^(^a 4“ ^oo) — 

^2 — ^(^a ^oo) — ^45 

2 (A.l) 

= “ 2(^0 + 6 * 1 ) = -hi, 

Vi = — 9 i) = — 62 . 

Note that the parameters 61 , 62 , 63,64 in (lA.lIl and the cr-function by Jimbo-Miwa, 
which we use all over the paper, are in fact those denoted in [46] as 6 ^, 6 ^, 6 ^, 6 ^ and 
6 +, respectively (see Eqs. (C.55)-(C.6I) of |4], and Eq. (1.13) and Props. 1.6 and 
1.8 of [46] 1. This means that the canonical variables q and p appearing below are 
related to those of P 6 given in Sect. 3.10 by a birational canonical transformation 
(see Eq. (1.12) of [46]). 


^The parameters of P6 used by Jimbo-Miwa in [4] and by Okamoto in m are related as 
9o = >co, Si = >ci, 9a = 9, and 0oo = J<oo + 1. 


38 





Given the set of parameters ,..., &4 and the cr-function, the canonical coor¬ 
dinate q and momentum p are given by the expressions 


9= ^[(^3 + ^4)5+ (o-'-6364)6'], 


1 


P = 


(cr' — W2)B -\- {wia — ui3)6]. 


(A. 2 ) 


2 Aq{q - 1 ) 

where 

A = (cr'-I-63)(cr'-I-64), B = £4(0; — l)cr"-I-Sicr'— S3, C = 2 (Q!Cr'— cr) — S2. (A. 3 ) 

Here, s^ and Wk, k = 1 , 2 , 3 , are the fc-th fundamental symmetric polynomials of 
{61,62,63,64} and {61,63, 64}, respectively. The corresponding system of Hamilton¬ 
ian equations reads 


9 = 
p' = 


9 ( 9 - l)(9-a) 


a{a — 1) 
1 


o ^0 di 
2p -- 

9 9-1 


q — a 


y{a — 1) 


I [- 3 ( 3 '^ -I- 2 (q; -I- l)g - a] p^ 


+ [6*0(29 - 1 - a) -I- 6 *i( 2 g - a) + 6 a{ 2 q - l)]p- 


(00 + 01 + 0a)^ ~ 


(A. 4 ) 


Note, that eliminating p, one finds that q solves the canonical form of P6 with 
a = |0^, 6 = —^00, c = ^ 9 i, and d = —^0a(0a + 2), in agreement with the 
birational canonical transformation (which can be identified as the change 0^ 1—>■ 
0 a -I- 1 , 0 OO 0 OO + 1 in the set of monodromy parameters of P6 given in Sect. 3 . 1 ). 

Now we are ready to consider our case. In Sect. 4.2 we impose the following 
relations on the parameters: 


r — s 


Vl=V‘i = 


V2 = —= 


2 ’ 2 
Hence, 63 = —61 = —ui and 64 = 62 = V2, and so we have 


Si = 2 z/2, S2 = -vi 

Wi = V2, W2 = -1^1, 

We find it convenient to set 


S 3 = —2^11/2 
W3 = -h'lv2, 


Vl = SUi, 


where 


1 


V2 = SV2, 

I — V 


Vl = -- 


2 v 


V2 = 


2 v 


We recall that s is a large parameter, while vi and V2 are finite. 

Let us map the system (IA. 4 I) into the form considered by Wasow. Consider the 
formal asymptotic expansion for the cr-function (I 4 . 15 |) . which is uniquely defined 
by its first term, 

■ +v'^ \ 

—-—- — 2 viV 2 \fa I -b 0(1), s ^ 00 . 
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Using it and equations (jA.2l) and (IA.3I) . one finds the corresponding formal expan¬ 
sions for q and p in 1/s. The first terms reads 


q = —va + s 


^ {a-l)y/a{vi+V2) 

4,{vl^/a- V2)(v2^/a- vi) 


+ 0(s-^) 


Viy/a-V 2 , (■v/a-l)[(wi 

p = s ^ -I- 




(A.5) 


^/a[^/a + l) 4^/a{y/a + l){v2y/a - vi) 
Introduce a 2-component vector function y = y{a), 


y = 



9 = 2/1, P = sy2, 


where this transformation is motivated by the formal expansion (IA.5I) . Taking into 
account that in our case do = —Goo = —vi + 1'2 and 9i = —9a = —vi — 1^2, can 
write the Hamiltonian system (IA.4I) in the following form: 


e 



f{oi;y), 


where e = 1/s and 


/(a; y) 


( h{a\y)\ 

\f2ia-y)) 


(A.6) 


fi{a',y) 

f2{a;y) 


2 / 1 ( 2 /! - 1 )( 2/1 - «) (a Vl - V 2 _ iyi +V2){0‘- 1 ) \ 
a[a-l) V yi [yi - l){yi - a)) ' 

—TT{ [ - %! + 2 (a -f 1)1/1 - a\y2- 2 (ui - 112)2/1 + 2 i;ia - 2112}. 
a{a — 1) ■' 


The Wasow theorem deals with the system of the form (IA.6ll^ and contains two 
additional conditions, called in |55j assumption (A) and assumption (B). We are 
going to verify that these conditions are valid in our case. 

To verify assumption (A), see p. 218 in |55) . we first assume that a € D where 
D is any closed disk inC\{0,/3,l,oo}. The case of our particular interest is where 
a belongs to any closed segment in the interval (/3,1). Assumption (A) requires 
also existence of a holomorphic vector function (f> on D such that 


/(a; $) = 0. 


In our case 




4>i = -y/a, (j)2 = 


Vly/a — V2 

v^(v^ + 1) 


(A.7) 


Note that is holomorphic for a £ D. Thus, assumption (A) is verified. 

To verify assumption (B), see p. 219 in [55], we define a new vector function 


y* = y 


and rewrite (IA.6I) as follows: 


d ^ d - - N 

e^ 2 /* = + -^y* + d(«; y*)- 


(A.8) 


Equation (IA.8I) exhibits the constant (with respect to 1 /*), linear and nonlinear 
terms of the system of equations (jA.6|) about the zero-order (in e) solution (/; the 


^Note, that here we meet a special case of the Wasow theorem, where f{a, y) is independent 


of e. 
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vector function g{a]yt:) satisfies dy^ y*)\v,=o = 0, fc = 1,2. Equation (IA. 8 I) is 
a special case of Eqn. (36.9) in [55]. The matrix A reads: 


A = 



a(l - i/a) 

[uiq:( 2 + i/a) - r)2(l + 2i/a)] {vi^/a- V 2 ) 
a'^(l — a)(l + y/a)'^ 


q I + \/« \ 
1 - ^/o^ 
Vly/a-V2 
a(l - ^/a)/ 


The determinant evaluates as 


det A = 4: 


{vi^/a- V2){v2y/a 
a3/2(l-a)2 


^’i) 


Recall that the disordered regime which we consider in Sect. 4.2 corresponds to 
a € (/3,1], where /3 = = v^/vl- Thus, det^ ^ 0 for all values of a of this 

regime, and assumption (B) of the Wasow theorem is fulfilled. 

Hence, the Wasow theorem applies to the system (IA. 6 I) , since both assumptions 
(A) and (B) of the Wasow theorem are satisfied in the disordered regime. Therefore, 
from the Wasow theorem it follows that there exists a genuine solution q, p of 
the system (|A.4p with the asymptotic expansions (jA.5|) . These expansions are 
uniform on any closed segment of (/3, 1 ) and all power-terms in e can be constructed 
uniquely via the explicitly given zero-order solution (IA.7[) . Due to the one-to- 
one correspondence between the canonical variables 9 , p and the cr-function, see 
Eqn. (1.6) and Prop. 1.2 in [46] . there exist genuine solution cr-function with the 
formal expansion (14.151) . 

Note that the switch to the Hamiltonian system q and p has mostly theoretical 
rather than practical sense since from the construction above it is not obvious that 
the resulting expansion for the cr-function contains only even powers of s. However, 
this fact is obvious due to the cr-form of P 6 (II.3p . 


Appendix B. Saddle-point approach to asymptotic expansion of the 

EFP in the ordered regime 


As in the main text we assume here that <7 = 0, which is not crucial for the 
approach discussed below. As mentioned in the main part of the paper the complete 
asymptotic expansion of Fr^s,q can be presented as transseries with a finite number 
(actually, of s) series. The transseries expansion of logFr-jS.g consists of infinite 
number of series; the n-th series in this expansion is generated by the term Tr(A"’). 
Here we consider only construction of the first series of the complete transseries 
expansion, i.e., the series generated by Tr(Ar); the other series can be studied 
analogously, however the corresponding calculations are more involved as they are 
related with the saddle-point method for integrals with n -I-1 variables. 

The main object of our studies here is the following integral: 


Tr(A:) 


a 


r—s+l 


(2^)2 



1 (1 — lyy {1 — pY du dp 

avpY (1 — avY (1 ~ apY v^p^ 


(B.l) 


We rewrite it in the following way 


Tr(A:) 


^r-s+l 


(2^)2 


I I e^SMf{^,p)di,dp, 

J Co Co 
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where 


and 


S{iy,fj.) = S{iy) +S{n), 


S{iy) = - log Y —— - log(-i^), 
V 1 — aiy 


(B.2) 




1 


(1 — au^Y ’ 

Note that in the definition of we used the fact that s is an integer and we 

recall that v = s/r. 

Since the variables v and fj, are separated in the function one can 

successively apply the usual saddle-point method for a single variable with respect 
to each of these two variables. However, it is easier to employ the saddle-point 
method for multiple integrals. We follow below the version of this method developed 
by Fedoryuk [HT] . 

Saddle points are defined by the system of equations -^S = -^S = 0. There are 
four saddle points with the coordinates (z/_,p._|_), (z/_|_,/r_), and 

where 

i^± = /i± = 

In the following we assume that 

0 <v <u= ]~^ (B.3) 

1 + Va 

This inequality, together with the condition 0 < a < 1, implies that 

<v- <Q. 

To prove it one has to notice that v+v- = 1/a where both factors in the left-hand 
side are negative. Using this relation we find that 



1 — V 

S{v-) + S{v+) = -logo > 0. 


(B.4) 


If we denote y = —v and consider ■^S{—y) for ?/ > 0, then it is easy to notice that 
-^S{-y) < 0, for 0 < 2 / < ?/_ and y > y+, and ^S{-y) > 0, for y_ < y < y+, 
where y± = —v±. Therefore, 

S{v+,y.+) > S{y+,y-) = S{y-,y,+) > S[y-,y.-). 

In other words, the saddle point surface should not contain the points (z/+,/i+), 
(i/_,2i+), and (v+,y._). 

Thus, asymptotic expansion consists of the contribution related with the sad¬ 
dle point We use the general formula for this contribution obtained by 

Fedoryuk m. which in our case can be written as follows. 


Q,r-s+l ^2sS{v-) 




(B.5) 


where 


Rk{s;a,v) = 


k\ 


( ^ ^ 

k 

\ ( 


1 1 


(/(u m)' 
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and A is the Laplacian, 


dv'^ d^jL^ 


The function /x; /i_) is 




S{i^, fi) - S{iy-,n-) - {{ly - v_f + - fi-f) , 


where 




(1 — v ) v - + 1 + u 
vlil-v.) 


It is important to mention that (IB.51) produces exponential decay for large s. To 
prove this, we rewrite identity (IB.4p as 


S{v-) + S{v+) + ^^ log a = 0 

V 

and employ the inequality S{v-) < S{v^), that yields 

2S{v-) + ^^ log a < 0. 

V 

It is interesting to note that S{v-) > 0. To show this let us consider S{—y) for 
positive y. Since for 0 < y < 1 both terms (see the second equation in (IB.2I) 1 are 
positive, S{—y) > 0 for these values of y. At the same time 

5'(-y) =-logy--logo+ 0(l/y) < 0, y-)>oo, 

V 

so that all positive zeros of S{—y) should be greater than 1. On the other hand, 
it is easy to see that S(l/y/a) > 0, therefore one zero is greater than ?/+ and the 
others (if any) should be located in the interval (1, To prove that actually 

S{—y) > 0 for y G (1, Xj%Ja) consider the following change of variables, y = l/a'^, 
where 0 < e < 1/2. In this notation S{—y) reads 


S{-y) 


-log 

V 


1 + a" 
1 + 


1 - r;, 1 „ 

-log — >0. 

V 


Now we study expansion (IB.5|) in more detail. Let us note that (IB.5|) does not 
represent an asymptotic series because of differentiation of the function exp(sS'i) 
with respect to the large parameter s. To get the asymptotic series one has to 
rearrange summation in a proper way, such that the series would represent an 
expansion in the decaying powers of s. Let us write, 


[f] 


Rk { s ; a , v ) = ^ 


1=0 


Rk,iia,v) 

gk-i ’ 


where [^] is the entire part of Note that 
Ro{s;a,v) = Ro,oia,v) = 
Thus the asymptotic expansion of Tr(Ar) reads 




1 — V 

V 




bp 
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where 


In particular, 



1 


2p 

'^Rp+i^i{a,v). 

1=0 


bo — 1 , 


h 

h 


Ri,o{a, v) + R 2 ,i{a, v) + Rz, 2 {a, v) 

R 2 ,o{a, v) + ^3,i(«, v) + j?4,2(a, + Rp^a, v) + ^6,4(«, v) 


To evaluate &i, we note that 


i?uo(Q;,w) 

R2,iia,v) 


■R3,2(a,w) 


12 S"{v-f 


where 



d 


and we use the convention that Df{v-,fj,-) = /i)_ 

above formulas, we obtain 


Using the 


(1 — a)'^{v‘^ + 1) + 9(1 — a‘^){v^ + v) — 2(10a^ — lla + 10)u^ 

6 vT ^-^((1 — vy — (1 + 


Similarly, to compute 62 , we note that 


R2,o{oi,v) 

Ro,i{a,v) 


Ri, 2 {a, v) 




Ro^4{a,v) 


8S''{i2_y ’ 

+ 765(3)45(3)(^_)2(i7A + 3D^)^f{iy-,fi-), 

(ll5W(._) + 85(3)(._)i^) /(._,^_), 

205 5(3)(^_)4j(^_^^_) 

288 5"(:/_)6 
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Using these formulas we obtain 


^2 


■ 72(l-a)[(l- 
+ 18(l + a)(l 


b - , ^2 13 [(1 “ + 1 ) 

v)^ — (1 + vYa\'^ 

af{v'^ + v) + (llSa^ + 782a + 113)(1 - af{v^ + v^) 
18(1 - a^)(35a^ - 36a + 35)(v^ + v^) 

+ 12(83(a'‘ + 1) - 194(a® + a) + 321a^)i;'‘]. 


This result can be also written as follows: 

bi — bi , ^2 ^ “t“ ^2 ■ 

Here hi and 62 are given by (I4.36|) . where they are rewritten in terms of u (see 
Equation (IB.3l) b 


Appendix C. Explicit form for the exponentially small correction in 

the ordered regime 

As explained in Appendix [B] the complete asymptotic expansion in the ordered 
regime can be written in terms of transseries. The first series is defined by Tr AT and 
given by the double integral (IB.ll) . At the same time our construction of asymp¬ 
totics presented in Section 14.31 makes it possible to obtain another representation 
for this series. We recall that in terms of the a-function the corresponding series 
is written in (14.281) as where the function w is a solution of ODE (14.301) . The 
general solution of (|4.30|1 can be found explicitly, namely, it has the following form 


SUJ = 


— b— {z"^ — d)a 
where the function z reads 


l-\- V 


, b = 


1 — u 


2 a d 

z = ± — —logip, 
s da 


with 


\1-I-s(l - u)/u ) 




(C.l) 


(C.2) 


s, —sjv 
— s(l — v)lv 


a 


Here Ci and C 2 are the constants of integration, which, in fact, are functions of s 
and V. 

Now we have to choose the solution which corresponds to Ei-_s,o- As follows 
from (|4.14l) and (j4.28|l . 

logE,,,.o = - Te-®-—^^ + 0(e-2®-). (C.3) 

Jo a(l -a) 


Using asymptotics as a 0 of the Gauss hypergeometric function and comparing it 
with the expansion (14.251) we find that C 2 = 0 and can fix the constant of integration 
while finding w from the first equation in (IC.lj) . Before presenting our final result 
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we give an intermediate formula for suj' which can be obtained by manipulation 
with (1C.II) and (|C.2|) at C 2 = 0, 


da 


log 


-h-{z^ -h)o^ 2 F 1 


s/v, —s 
1 + s{l — v)/v 


After integration and substitution of sui into (1C.31) with yet indeterminate constant 
of integration, application of the formulas for differentiation of the hypergeometric 
function [56], determination of the integration constant as explained above, and 
some further simplifications, we arrive at the following representation of logfr.s.o 
in the ordered regime. 


logFr,s,o = 

s(l — v) 


1 — a 


s(l — ■!;) + ' 


■ 2F1 


s/v, —s 
s(l — v)/v 


2FI 

a I 2 F 1 


s/v, —s 
1 + s(l — v)/v 

1 s/v, 1 — s 

2 + s(l — v)/v 


n 2 


da 

+ O (e-2®“) . (C.4) 


To see that this integral, in fact, is exponentially decaying as s —>■ 00 , one has 
to apply the saddle-point method to the corresponding integral representation of 
the hypergeometric function. Contrary to considerations of Appendix jB] the as¬ 
ymptotic analysis of the representation (IC.4|) requires just the saddle-point method 
applied to a single integral. Note that the correction term in (1C.41) has the order of 
the square of the first term. 
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